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Abstract: Non-extremal solution with warped resolved-deformed conifold back- 
O , ground is important to study the infrared limit of large N thermal QCD. Earlier 

works in this direction have not taken into account all the back-reactions on the 
geometry, namely from the branes, fluxes, and black-hole carefully. In the present 
^ . work we make some progress in this direction by solving explicitly the supergrav- 

ity equations of motions in the presence of the backreaction from the black-hole. 
The backreactions from the branes and the fluxes on the other hand and to the or- 
der that we study, are comparatively suppressed. Our analysis reveal, among other 
things, how the resolution parameter would depend on the horizon radius and how 
the RG flows of the coupling constants should be understood in these scenarios, in- 
cluding their effects on the background three-form fluxes. We also study the effect of 
switching on a chemical potential in the background and, in a particularly simplified 
scenario, compute the actual value of the chemical potential for our case. 
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1. Introduction 

Much of what is known about the phases of strongly coupled gauge theories (and 
in particular, QCD) comes from a variety of techniques, each of which accompanied 
by its attendant limitations. Perturbative (i.e. weak coupling) computations can 
probe a large part of the parameter space of the theory, like allowing one to deal 
with varying number of colors iV, flavors Nf. However, these results are valid only at 
temperatures well above the deconflnement temperature Tc, and at large values of the 
baryon number chemical potential /i in order for the QCD coupling to be small, and 
thus the perturbation valid. These exclusions put almost all of the interesting region 
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of the parameter space explored by RHIC data beyond the reach of perturbative 
computations. 

Lattice gauge theory, which provides a rigorous non-perturbative starting point 
for QCD, is not without its hmitations as well. It is difficult to incorporate realis- 
tic quark masses, and results are limited to the regime near Tc, and /i very small. 
Nonetheless, a combination of such conventional methods of analysis (including in- 
sights from effective theories like chiral models) suggest that the gauge theory pos- 
sesses a color superconductivity phase at asymptotically large value of the baryon 
number chemical potential /x. The literature is replete with conjectures for the phase 
diagram of QCD in the (T, /i) plane, especially for small values of T and large values 
of yU (see for example 

In recent years, there has been a considerable advance in understanding the be- 
havior of U{N) gauge theories at finite temperature using the gauge/gravity duality. 
That this development is more than timely is beyond dispute, as the new and inter- 
esting results from RHIC have provided a glimpse into a wide variety of interesting 
phenomena arising in the strong coupling regime of QCD. For instance, the quark- 
gluon plasma (identified as a new state of matter) displays many properties of a 
fiuid with low (shear) viscosity, explanations for which are difficult to obtain from 
traditionally available tools in perturbative QCD. 

Many (perhaps most!) analytic results coming from gauge-gravity duality are 
derived for gauge theories with TV = 4 supersymmetry and with N very large, and 
in the limit of the theory possessing exact conformal symmetry. One may thus 
genuinely be concerned about their applicability to QCD for which all these are not 
true. Recent progress in this area, however, has provided us with strong hints to 
overcome these limitations, and move towards models of gauge-gravity duality that 
are not supersymmetric, and are non-conformal (in a sense that will be made precise 
later). 

The first model that managed to expand the original AdS / CFT construction to 
incorporate renormalization group runnings is the Klebanov-Strassler (KS) model 
(with an extension by Ouyang ||25| to incorporate fundamental matters) that provided 
an IR dual of, although not exactly QCD, but at least its closest cousin: large 
supersymmetric QCD. The UV of the original Klebanov-Strassler model is now 
known to have some issues, like the divergences of the Wilson loops at high energies, 
and additional Landau poles once fundamental matters have been introduced. This 
means that UV completion is necessary, and to have the full gravity dual of the 
corresponding gauge theory that behaves well at high energies, the KS geometry has 
to be augmented by a proper asymptotic manifold. 

Other extensions to the original KS model quickly followed. For example in 
the cascading picture of the original KS model was extended to incorporate black-hole 
without any fundamental matter, which was then further extended to incorporate 
matter in [H]. However none of the above models actually considered the full UV 
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completion as most of the analysis of these works were directed towards unravelling 
the IR physics. Therefore issues like UV divergences of Wilson loops and Landau 
poles were not investigated. 



In a series of works 16, 20, IS] done over the last couple of years, we tried 



to address these concerns. Our aim therein was to incorporate the backreactions 
from the black-hole, fluxes, and branes consistently so as to have a well defined UV 
completion that not only allow us to get rid of all the poles etc., but also give us 
a model that could come closest to what we might have expected from a large 
thermal QCD. We did manage to at least successfully generate such a UV completed 
dual picture, but many of the backreactions turned out to be too difficult to incor- 
porate fully. One aim of this paper is to make progress in this direction. In sec. 



2.1 we will show how exactly to incorporate the backreactions from the black- hole 



in IR regime of our theory to lowest orders in string coupling, and color-to-flavor 
ratio. Interestingly, to this order, the backreactions from branes and fluxes could be 
consistenly ignored. We will demonstrate this in sec. |2.1.1| , and the full EOMs will 
be solved in sec. [2.1.2| and in Appendix |C[ 

Although to the order that we study the IR regime of our theory allows us to 
ignore the backreactions of the fluxes, we will in-fact work out the detailed fluxes in 
sec. E.21 The backreactions from the black hole and the flavor seven-branes will be 



fully incorporated in the fluxes. It will also be clear from sec. 2]2| as to how these 
backreactions conspire to make the three-form fluxes non-ISD. 

One persistent problem associated with thermal QCD is the interpretation of 
beta function of the theory. There is a long history on the subject starting with the 
classic work by Collins and Perry |]TU|. In sec. ^.3| we discuss briefly how we should 
interpret the running of couplings from the gravity dual. Our work strongly suggest 
one framework, although alternative interpretation could possibly be entertained. 

We end sec. ^by taking a short detour. In sec. ^]4|we study the dipole and non- 
commutative deformations of the seven-brane theory on the gravity side. This detour 
is not without its merit. Dipole deformations are known to increase the masses of 
the KK states, so that a dipole-deformed theory will have a slightly different KK 
spectra. We again make very brief speculations of the underlying physics, leaving 
most of the details for future work. 

Thus sec. prepares us with a backreacted metric, and with backreacted fluxes. 
One may now go on with this to flU up rest of the missing steps left in R |2D|, ff^ 



These issues however will be addressed in future works. Here we aim for more modest 
return. In sec. |^ we present a small computation on the chemical potential. The 
reason for choosing this computation over other possible interesting ones is two-fold. 
One: its simple enough computation that carries sufficiently interesting physics, and 
two: we use this to show in sec. |3.2| how an alternative way of getting the chemical 
potential, via say duality chasing, may be inherently flawed. In sec. p.3| we compute 
the chemical potential for our model. In Appendices |^, and ^ we speculate more 
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on the duality chasing techniques. We conclude in sec. ^ with some discussions. 



2. Analysis of the background 

Let us start with the model that we studied in 0. The IR physics is captured by 
the Ouyang-Klebanov-Strassler-black-hole (OKS-BH) geometry, namely, the small 
r physics is determined by a warped resolved-deformed conifold with fluxes, seven- 
branes and a black hole in the ten-dimensional spacetime. On the other hand the 
UV physics is conformal, and is captured by an asymptotically AdS geometry with 
fluxes and seven-branes. 

As discussed in [|l6l, these two geometries, namely the asymptotic AdS and OKS, 



can be connected by an intermediate configuration with brane sources and fluxes. 



These branes sources were elaborated in details in [16|, although many coefficients 
in the background geometry were left undetermined therein. In the following we will 
fill up some of these missing steps. 

Let us begin with the basic ansatze for the metric in the three regions. For all 
the three regions we assume that the radial coordinate r spans b < r < rmin for 
Region 1 where we expect all the confining dynamics to take place; rmin < r < Tq 
for the intermediate region called Region 2; and < r < oo for Region 3 which 
captures the asymptotically conformal region. The minimum radius r = b, which 
signifies the cut-off coming from the blown-up (as well as S^, although for most 
of the calculations in this paper we will only consider a warped resolved conifold 
instead of a warped resolved-deformed conifold), maps to the expectation of the 
gluino condensates of the dual gauge theory at zero temperature. Considering all 
these regions, the non-extremal metric takes the following form: 



1 



gi{r)dt^ + dx^ + dy'^ + dz^ + Vh g2{r)'^dr'^ + dMl 



ds = —= 

= _e2^+2B^^2 ^ e^A^.^dx'dx^ + e-'^'^-^^mndx'^dx'' (2.1) 

where gi{r) are the black-hole factors and we have taken gi = g2, the components go 
as i,j = 1, 2, 3 and m,n = 4, 9, the warp factors A, B are defined as: 

A = ~ ^log B = hog g, (2.2) 

dAif is typically the metric of warped resolved-deformed conifold and h is the warp 
factor that behaves differently in the three regions as shown in [|l^ . 



Observe that in the extremal limit, (^i = 5'2 ~ 1 and the extremal metric is dual 
to the low temperature confining phase of the gauge theory. To see this, note that 
in the absence of any seven branes. Region 1 of the geometry of |16[ in the extremal 
limit is identical to the IR geometry of Klebanov-Strassler (KS) model 0. If seven 
branes are placed far away from Region 1, that is ^ b, we can neglect their 
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back-reactions and consider the axion-dilaton field to be effectively constant as in 



p. Hence in the extremal limit, Region 1 of is identical to the IR region of KS 
which, in turn, is dual to the low temperature confining phase of the SU{M) gauge 
theory wherein chiral symmetry is broken. The extremal geometry can incorporate 
temperature of the field theory once we analytically continue to Euclidean signature 
with it ^ T and impose periodic and anti-periodic boundary conditions for the 
bosons and fermions on the closed time circle. Furthermore, in extremal case the 
entropy will vanish. This is expected as the entropy from the dual geometry arises 
from the fluxes which are at least 0{Nes), where N^s is effective brane charge. As 
the deformed cone represents confinement of charge, we expect to get A'^efr = from 
the dual geometry. This is indeed what happens as energy scale for a thermal field 
theory is set by the temperature and at low temperature, only the IR degrees of 
freedom are excited. This means in the dual geometry, all we need is the region near 
r ~ 6 of the deformed cone — but in this region the five-form flux vanishes and 
we get iVeff = 0. 

As the temperature is increased, we expect that the non-extremal solution will 
have less free energy than the extremal solution, just as in the case for the AdS-black 



holes 10], and Hawking-Page phase transition will take place [^. The focus of this 
work will be to analyze the non-extremal solution which is dual to the deconfined 
phase of large N thermal QCD, while a detailed analysis of phase transitions will be 
presented in a follow up paper ||12 . 



The non-extremal solutions we present in this paper are precisely dual to the 
high temperature regime of the gauge theory — where chiral symmetry is restored 
and the light degrees of freedom are deconfined. However, heavy quarkonium states 
arising from the seven branes placed in the UV region can coexist with the chirally 
sjTumetric phase above the deconfinement temperature. But as temperature is raised 



even further, the heavy quarkonium states will eventually melt |20|. 

For both extremal and non-extremal cases, typically h would have logarithmic 
factors in Region 1 whereas it would have inverse r behavior in Region 3. In the 
intermediate region, the warp factor will typically have both the logarithmic and 
the inverse r behavior. Therefore to summarise, the background should satisfy the 
following properties: 

• Fluxes are non imaginary self-dual i.e non-ISD, and become ISD once the black- 
hole factors Qi in the metric are removed. Therefore the deviation for ISD property 
is proportional to the horizon radius Vh- 

• The gravity dual of the deconfined phase is given by resolved warped-deformed 
conifold with a black-hole. In the limit where the deformation parameter is small, 
the background can be succinctly captured by a resolved conifold with fluxes and 
black hole. 
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• The resolution parameter is no longer constant because of the various back-reactions. 
In fact the resolution parameter becomes function of r/j/r as well as QsNf, and 
QsM^/N where Qs is the string coupling, M is the number of bi-fundamental matter, 
is the number of colors, and Nf is the number of fundamental flavors. 

From the above set of arguments, we can use the following ansatze for the internal 
metric: 



QfYiff^dx doc 



dr + r e 



2 2B 



A^dip + cos 9id(f)i + cos 62 



+ 0{gsM^/Nyjr') 



B{ddl + sin^^i^^) + (1 + F){1 + g) 



6 



+ 2/b cos ip{d9id92 + sin 61 sin 62 



"^^2 +sin2^2rf02 



(2.3) 



hd4'2) ~ sin iplsin 62d(f)2d0i — sin 6id(j)id62 



where we will only consider the resolved conifold limit, with F being related to the 
resolution parameter (whose value will be determined later). In other words, we take: 



/, ^ 0, F ^ g ^ 



A 



B 



o 



(2.4) 



where the numerical factor of 6 is inserted to bring certain expressions in a better 
format. As we will see, this F (or equivalently a) determines the squashing factor 
between the two spheres, and we can consistently keep the second squashing factor, 
to be zero. 

In the following section we will argue for these parameters and their dependences 
on the horizon radius by analysing the non-extremal limit of the warped resolved- 
deformed conifold background^. 



2.1 Derivation of the non-extremal BH solution for the Klebanov-Strassler 
model 

We first compute the non-extremal metric arising from Type IIB supergravity action 
given, in the notations of [r3|, in the following way^: 



1 



IIB 



+ 



2^10 



R- 



2|ImrP 12Imr 4 ■ 5! 



C4 A G3 A G3 
Imr 



5i 



loc 



(2.5) 



^We will continue calling this background as the Klebanov-Strassler background as they all fall 
in the same class of supergravity solution. 

^Although in this section we will use the Einstein frame to express the metric, we will however 
not distinguish between the two frames in later sections because the dilaton will be considered 
constant, unless mentioned otherwise. 
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where Sioc is the action for all the localized sources in ten dimesional geometry 
i.e five-branes and seven-branes mostly from Region 2 onwards. Our aim is to re- 
analyse the non-extremal Klebanov-Strassler solution. Recall that for Klebanov- 
Tseytlin model the non-extremal solutions were analyzed in [1^1 , while in ^ there 
have been studies of gravity duals of finite temperature cascading gauge theory with 
fundamental matters^. However in |^ precise background fluxes and the warp factors 
taking into the backreactions of the BH geometry were only conjectured. Here we 
will derive the non-extremal metric dual to a UV complete gauge theory that mimics 
features of large N QCD at the lowest energies, justifying the proposals made in 
One immediate outcome of this would be the verification of the conjectured 
dependence of the resolution parameter on the horizon radius r^. 

Our ansatz for the metric is ( |2.1| ). We look for solutions with regular Schwarchild 
horizon at r = Vh- This is achieved by imposing e^^^'*-' = and considering solutions 
to A such that e'^^^'^^ ^ 0, which guarantees a non-sigular horizon By solving 
the Einstein equations along with the flux equations with these boundary conditions, 
we will find the non-extremal solutions with regular horizons. 

Observe that we have warped Minkowski four directions, a non-compact radial 
direction r and a compact five manifold M.^. The back reactions of the fluxes Gs, 
and axion-dilaton field r will modify the warp factor e"^"*"^ while 'gmn will be al- 
tered due to the presence of a black hole and the various sources. In particular gf^n 
will be a warped resolved-deformed conifold with a bare resolution parameter a^. 
Note however that only the warp factor e^^^ will be essential to analyze the con- 
finement /deconfinement mechanism for the boundary field theory |T^. The linear 
confinement of quarks and the string breaking mechanism which eventually describes 
the deconfinement of QQ pair, is only sensitive to the warp factor. The exact solu- 
tions for the internal metric in the non-extremal limit taking into account the back 
reaction of the various fiuxes is not essential to study free energy of the QQ pair. 
Nevertheless we will find the exact form of the internal metric upto linear order in 
resolution function F . 

We restrict to fluxes and axion-dilaton fleld r which only depend on and not 
on the Minkowski coordinates . Then the Einstein equations can be written as 



R 



R-n 



G3 ■ G3 
48 Imr 

G3 ■ G3 
48 Imr 

mn o y 



+ 



+ 



It 

8-5! 



It 
8-5! 

rploc 



+ 



+ 



/labcd-' 



' abed 



4-4! 

rp rp abed 

^ mabcd^ n 



+ 



T^loc rp\{ 



4-4! 



+ 



4 Imr 



nbc djYiTOnT 



2 Imr 



(2.6) 



^See also [nsi where somewhat similar analyses were also done. 
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where is given by the following self dual form 

F5 = (1 + *io)da A A dx^ A dx^ A rfx^ 
with a = e^^ and T^""^ being the trace of 

2 6 Sloe 



loc 
ab 



-9 



ab 



Using the form of the five- form flux (|2.7| ), the first equation in (|2.6| ) becomes 



R 



G3- G3 



loc 



(2.7) 



(2.8) 



(2.9) 



48 Imr 4 

On the other hand, the Ricci tensor in the Minkowski direction takes the following 
simple form 

1 



R 



dmig'^'^dng^u) + ^?'""r*l,a„^7^, - g'^''g'''^'dmg>.'>.dngu'u (2.10) 



where = 0, .., 3 and F*^^ is the Christoffel symbol. Now using the ansatz (|2.1| ) 
for the metric, ( p.lOj ) can be written as 



Rtt = e^(^+^^ [y\A + B)- 35^"9„59„(A + B 



where we have defined the Laplacian as: 



(2.11) 



(2.12) 



The set of equations can be simplified by taking the trace of the first equation in 
( [^.6[ ) and using ( |2.11D . Doing this we get 

V2(4A + B)- Sr^'dnBdU^A + B) = e-^A-2B GmnpG"''''' ^ e-'OA-A^'dmad^'a 

12Imr 



On the other hand using (2^) in ( p.ll ), one gets 

Ri -r: = q 

which in turn would immediately imply 

V^S - Sg^^'dmBdnB = 



loc 



(2.13) 



(2.14) 



(2.15) 

Minimizing the action ( p.5|) also gives the Bianchi identity for the five-form flux, 
namely 



dF^ = i/3 A F3 + 2k2qT3P3 



(2.16) 
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where ps is the D3 charge density from the locahzed sources [13|. Using ( p.7|) in 
( p.l6|) and subtracting it from ( p.l3D one gets the following 

V2(e4^+^ -a) = ^KG3 - H Gs\' + e-^^-^^|a(e^^+^ - a)\' 
olmr 

The Ricci tensor, on the other hand, for the x"^,m = 4, ..,9 directions takes the 
following form 

Rmn = Rmn + Qmn^^ + B) - Zg^nQ^^ d^Bdu {A + B) 

+ 3V„9„5 + dmBdnB - %dmAdnA - 29(^A9„)5 (2.18) 

where is the covariant derivative given by 

V^K = a^K-fLH (2.19) 

for any vector V^,. Here Rmn is the Ricci tensor and V^^^ is the Christoffel symbol for 
the metric 'gmn- The equation for Rmn is given by: 

p _ _ G^- Gmab ■ G'^ drnXdnT 

^mn Qmn r, ^ t ~r TT ~r 



24Imr 4Imr 2 | Imr ^ 

P rpabcd fp rpfiabcd 
H 7j H fifmn 77 Tj H »Cm/lCrt/l 

4-4! lo ■ 4! 

-3V^9„S - 9„59„S + 29(^A9„)5 (2.20) 

which means, in general, this could lead to twenty different equations in six-dimensions 
(including another one for the trace). On the other hand the equation of motion for 
Gs can be expressed in terms of a seven-form A7 = *ioG3 — iC^ A G3 in the following 
way: 

dkj + A ReA7 = (2.21) 

Imr 

where typically A7 would study the deviations from the ISD behavior. For example, 
using our metric ansatz we can express A7 as 

A7 = [e^^+^ *6 - mGs] A dt A dx A dy A dz (2.22) 

The above choice of A7 leads us to three different classes of solutions from the G3 
EOM ( p.21|) . These three classes can be tabulated in the following way: 

• If a = e^^+^ in ( ^1221) and A7 = rfA7 = then G3 must be ISD. When B = then 
this is the same as GKP solution [T^, and in this case r is not restricted. 

• If a 7^ g4A+B ^i^QYi we can take A7 7^ but keep dAj = and dr = 0. This means 
A7 = c/Afi, where Ar is a six-form, and r is a constant^. 



^Or T = dX_i i.e d of a (— l)-forin. The functional form for the (— l)-forni is non-trivial, so this 
option is more cumbersome to use. 



- 9 - 



• If a 7^ g4A+B j^i^Qj^ (,gj^ again take A7 7^ but now dAi ^ and dr ^ such that 
( p.21|) is satisfied. This means both axion and the dilaton could run in this scenario. 

In this paper we are taking a = e^"^, so we have to consider the last two cases. 
Expressing A7 as Aj = A dt A dx A dy A dz we have *6 G3 — iG^ = T3 where 
T3 is non-zero as long as B is non-zero. A simple solution then would be to restrict 
oneself to the second case, i.e 

dT^ = 0, r = constant (2.23) 

Notice also that at far infinity, i.e r — )• 00, i? — )■ 0, therefore T3 — )■ as well^. Using 
the above argument, G3 can then be expressed in terms of T3 as 

G3^q5il±^Ii (2.24) 

Since r = constant, this means the closure of will involve a non-trivial constraint 
connecting the internal metric components with B and T3. However in this paper 
we will not be solving these equations explicitly but approximating G3 by Ouyang- 
Klebanov-Strassler fiux Gf^ which is ISD in their metric. This approximation suffices 
for our case, as we show below. 

Let us substitute G3 = G^^^ and a = e^^ into ( |2.22|) . This will convert Ay to a 
simpler seven-form in the following way: 

A7 = e'^^(e^ *6 Gl - iG%) AdtAdx Ady Adz 

« 3e^^(e2^ - \)g''e^ahcGfde ^ d^'' ^ dx^ AdtAdx Ady Adz 

~ ^9""^rabfG^^l dx" A dx^ A dx^ AdtAdx Ady Adz (2.25) 

At large the right hand side is small and therefore deviation from OKS fiux 
is of 0{rh,l/N) so one may consider A « 0. This means G3 = G^^^ is a good 
approximation. Additionally, since F5 is self-dual, Rmn can be simplified as 



~ G3 • (jr3 Gmnh ' G'^ OmTdr,T . oTi-, ,-. . n . 

= -Srnn^ + + + 8(1 - e-^^)d^Ad^A 



-3V„9„5 - dmBdnB + 2di^Adn)B (2.26) 

We see the first two terms are suppressed by QgM'^/N and the third term is removed 
because r is a constant. So we can ignore these contributions for the time being. 
Then, assuming A and B only depends on r, (|2:26|) will lead to 



Err = 8(1 - e-^^)drAdrA - SVrdrB - OrBdrB + 2%v4a,.)5 



^This is of course without considering the UV completion. With UV completion the large r 
behavior is non-trivial as discussed in B M . 
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Rab = —drQabdrB (2.27) 

where (a, h) denote the angular directions. We now see that for r > r/i, the Rab 
contribution is suppressed equivalently as the Rrr contribution, therefore we need to 



keep both the parts. This conclusion can also be extended to Rmn in ( 2.18 ), which 
implies that the R^ and Rab contributions are equally suppressed. All this then 
further implies that we need to solve the twenty-one metric equations. This is a 
formidable exercise. Is there a way by which we can avoid doing this? 

A possible way out would be to study the relative suppressions of various terms 
in the system of equations. This criteria was already anticipated in For example, 
as we discussed in we can equivalently take: 

{gs.N,M,Nf) ^ (e^e-^e-^e-'^) (2.28) 

This would clearly show that {QsN, QsM) are very large but {gsNf, ggM^/N, g^MNf) 
as well as M/N are suppressed in the following way: 

{gsN,g,M) ^ (e^-^ 6^") 

{gsNf,gsMyN,glMNf,M/N) ^^c-d ^ ^c-2b+a ^ ^2c-b-d ^ ^a-b^ ^2.29) 
provided {a,b,c,d) satisfy the following inequalities^: 

a > b > c > d, a + c > 2b, 2c > b + d (2.30) 

Let the smallest scale in our problem be the ratio M/N. Then if the argument of the 
relative suppressions of various terms in Rmn has to make sense, one would require 
the precise range of r where our approximations hold water. This gives us: 



' ^ "-^{m) (2.31) 



Thus if we are in this range, we can see that the curvature terms simplify drastically. 
This would give us a hint that if we solve the simplest trace equation alongwith the 
flux equations ( |2.15| ), ( |2.16| ), and ( |2.17| ) we would be reasonably close to the correct 



answer because the other twenty component equations would only change the results'^ 
to (9(r^/r'^). So once we are in the range ( |2.31D the only corrections to our simplified 
trace equation will be to 0{gsM'^ /N) and 0{r\/r'^). This is not so bad because if 
we choose e in ( p.28| ) to be e = 0.1, then 



N = 10^ M = 10^ Nf = 10, g, = 0.0032, r > 17.78r/, (2.32) 



solution to the inequalities is a = 8, 6 = 3, c = 5/2, d = 1, as given in One can of course 
allow other values of (a, b, c, d) that satisfy the inequalities. 

^This in particular means that not only the coefficients of all the terms of the internal metric 
will change to 0{r^/r'^) but also any new component will appear to ©(r^/r"*). This is exactly how 



we choose our initial metric ansatze (2.3) and therefore the system is self-consistent. 
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which means for r beyond 17.78r/i the contributions coming from the individual 
component equations to the solution generated using only the trace equation will 
not be too drastic. 

Therefore, once the dust settles, tracing the second equation in (|2.6|), using 
( p7[3D , ( p7[5|) and (^), we get 



~mn 



DA 'XT 

b 6 o 



12|Imr|2 



= (2-33) 



where i? = g^^Rmn and we have ignored all local sources. 

Our goal now is to solve the system of four equations ( p.l5|) , (|2.16|) , ( p.l7| ) and 



( p.33|) and find solutions for the warp factors A,B, the internal metric Qmn and the 



fluxes. In obtaining the solutions, we will be working in the limit where there is no 
local sources, is closed while the explicit form of the fluxes that solve the flux 
equations are described in the following subsection^. As we mentioned earlier, if we 
choose a = e^^, {\i-l7\) will imply that is ISD, in the extremal limit i.e = 1. 



On the other hand, G^ is not ISD on a deformed cone in the presence of a black 
hole, and the terms in G3 which make it non-ISD are precisely proportional to the 
blackhole horizon and the deformation function F that appears in 'gmn- With these 
considerations and our choice of internal metric Qmn we get 



tGs-e^HGs\' = iTsp ~ 0(F^r«/r«) (2.34) 



Thus with a choice of a = e^^ + 0{F^), ( ^TtD can be solved exactly. But if F ^ 1, 



we can ignore terms which means upto linear order in F, ( |2.17|) becomes 

(2.35) 

Thus ignoring 0{F'^) in (|2.17|) ^, we are essentialy solving (|2.15|) , ( p.l6|) , (|2.33|) and 
( p.35|) . In fact we will show that (|2.16|) dictates F <^ 1 and our explicit numerical 



solutions will also be consistent with this assumption, justifying our perturbative 
analysis. 

Now only considering upto linear order terms in F, we get a = e^^ which relates 
the warp factor to the five-form field strength which in turn depends on G3 by the 



®It is of course possible to consider additional sources to obtain a UV complete solution as done 
in . But for the purpose of the current section, which is to analyze the non-extremal limit for 
the IR geometry, we will ignore local sources and discuss their effects briefly towards the end. 



^The term in (|2.34|) appearing in (|2.17] ) contributes as ~ 0{F'^{gsM'^ /NY),1 > 1 which can 
be easily obtained by using e"'''^ 0{gsN) [l + Oigs^P /N)]. As gsAP/N < 1, we can ignore 
0{F^gsM'^/N) terms. See also sec. |Tl]for more details on the various scaling limits. 
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Bianchi identity ( 2.16 ). Thus depends on the non-ISD G3 as G3 is modified in the 
presence of a black hole. But the choice of a = e^^ also means that the dependence of 
G3 on blackhole horizon rh appears in the form of a resolution parameter a = a{rh), 
a crucial fact that was first conjectured in ||^ and will be further illustrated in the 
next subsection. 

As already mentioned, equation ( |2.16D determining also depends on the in- 



ternal metric gmn- In the absence of any fiux and axion-dilaton field, gmn is the 
metric of base of the deformed conifold T^'^ which has the topology of 5*^ x 5*^. In 
the presence of a black hole horizon and various sources, the internal metric will be 
modified in the following way: 

9mn = 9mn ~^ 9mn (2.36) 

where Qmn is the metric of deformed cone with base T^'^ and therefore gmn denotes 
all the corrections due the black hole and all other sources. This means that gmn 
contains all the informations of the resolution factor and its subsequent dependence 
on the horizon radius etc. Note also that, as we have a horizon atr = rh with M units 
of fiuxes^° and Nf number of seven branes, gmn must at least be of 0{M, Nf,rh). 
We will evaluate e"^ and 'gmn to lowest order in and gsNf which in turn will 
drastically simplify our analysis. Our choice of gmn and gmn will be such that we 
have ( |2.3|) for the internal metric. 

The Bianchi identity for the five-form fiux, in the absence of any three-brane 
sources, reads 

dF, = H,AF3 (2.37) 

where F3 and are the RR and the NS three-form fiuxes. They are given as 

F3 = fP + 0{F), H3 = Hf^ + O^F) (2.38) 

where F^^H^ are the fiuxes in the absence of any squashing, that is for F = 
(we expressed this earlier as Cf^ = Fg^^ — rifg^^). For the regular cone, taking into 
account the running of the r field, F^ and are exactly the Ouyang fiuxes [EB], 
while the exact form of the fiuxes in a deformed conifold were discussed in 
Now from the form of the fiuxes on deformed cone^^ one gets that 

F3 ~ M[1 + C(F)], ~ gsM[l + OiF)] (2.39) 

Using this and ( |2.38| ) one readily gets that 

ifg A F3 = ^3^°^ A iff ^ + 0{M'^F) (2.40) 



^'^In the intermediate region, i.e Region 2 of the geometry, we will also have (p, q) five-brane 
sources. 

^^See section |2.2| for more details. 
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An immediate question is: what can be said about the squashing function Fl In the 
absence of the three-form fluxes, i.e M = 0, there is no squashing as the Klebanov- 
Witten solution with running dilaton needs no squashing. This remains true 



even when we introduce temperature. To see this, observe that the non-extremal limit 
of Klebanov-Witten(KW) model does not require any modification of the internal 
space: which means F = with e^^ = l — f^/r^ and the internal space is exactly T^'^. 
There could be squashing due to the running of r field in the KW blackground, but 
squashing would be at 0{glN'j), so we can ignore it as we will only consider up to 
linear order in QgNf. These behaviors indicate that F must be at least proportional 
to M. In the following subsection, we will justify this claim. 

2.1.1 Behavior of F and various scaling limits 

Let us go to the case when there is no blackhole but we have non-zero three-form flux 
i.e M 7^ 0. For this case we are back to Klebanov-Strassler-Ouyang background with 
no squahing and F = 0. This means, F must also be proportional to the blackhole 
horizon Vh- Combining this with the form of the Ouyang fluxes, taking into account 
of the back reactions of the seven branes, we expect 

F ~ 0{al, mgsM^lN^, r^^glN]) (2.41) 

with ao being the bare resolution parameter discussed earlier and (a, (3) are some 
integers. Notice that we have inserted a suppression factor of assuming /3 > 
in anticipation of a possible perturbative expansion. Therefore using our ansatz 
in (|O0| ) gives us 



F,^H, = fP a Hi'^ + 0{al rlglM-^'/N^^\ rlglN)) (2.42) 

implying that up to quadratic order in M, we only need Ouyang fluxes to solve 
( p.37|) . But to guarantee that we only need to consider upto quadratic order in M, 
we must show that higher order i.e 0{glM'^^'^ / N^^^) terms are small compared to 
the gsM"^ /N terms coming from Fg^^ A Hf\ This will indeed be the case once we 
solve (g]33) upto C(M2) We will see F ~ M/N where > M and this justifies 
ignoring the second term in ( |2.42| ) . In fact solving ( |2.37| ) with our ansatz for the warp 
factor shows that -^{F^ A H^) is the relevant term that enters into the equaton of 
motion (see Appendix 0). Hence in solving ( p.37| ) with our choice of warp factor, 
we are really ignoring 0{gsM^ /N"^) and keeping terms only upto 0{gsM'^ /N). This 
truncation is consistent for N ^ M which is achievable as we showed in ( p. 281) and 
( p.29|) . However one might question the suppression terms in ( p. 41]) and in ( p.42| ) 



12t 



If the solution to ( 2.37 ) upto 0{M'^) tells us that -F > 1, then we cannot ignore the second 



term in ( 2 .421 ) and therefore have to include 0{M^) and higher in solving (2.37). But, as we will 



argue soon, our solutions show that F 1, which justifies our truncation. 
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if {a,P) exponents are arbitrary compared to the range that ( 2.28|) would impose. 
That this will not be the case will become apparent from the following discussions. 
To start then we shall continue using the following five-form flux: 

F5 = (1 + *io)da A d^x (2.43) 

With this form of F5 and a = e^^ = 1/h, ( p.37| ) becomes an equation involving h, 
e^^ and F. We already know that in the AdS limit e^^ = 1 — f^/r^- In our non-AdS 
geometry we expect: 

^4 



1-4 + G 



(2.44) 



where G is at least 0{M, Nj). Using this expansion for e , alongwith the precise 
form of the Ouyang three-form fluxes F^^\h^^^ and only considering upto 0{M'^) 



terms ( p.37| ) reads 



(2.45) 



where (7™" is proportional to the deformed conifold metric (see Appendix jC]), h 
+ with being the Ouyang warp factor 




logr 



, SgsNf / 1 

ZTT \ Z 



H logr log sm — sm — 

M^N ^ ^ V 2 2 



(2.46) 



and is the contribution due to the presence of the black hole. 

We can readily see from ( ^^.45] ) why F ~ M/N . First note that the non-extremal 
limit of Klebanov-Witten model has an exact solution, h = L'^/r'^ with = 0. 
is only non-trivial due to the presence of three form fluxes, the black hole and other 
sources. Thus, ~ 0{M, gsNf,rh). On the other hand = QgNa'"^ and thus one 
gets from { ^J^ that F ~ 0{M/N, gsM'^ /N). But L^/a'^ > 1 and we can choose 
it large enough such that N ^ M which guarantees that F <^ 1. This is of course 
consistent with (OHI)^^. 



Again in ignoring higher order terms in M, we are assuming that F,G ^ 0{M) < 1, which will 
be consistent with our solution. On the other hand, the 0{M'^) term that enters into ( 2.45 ) from the 



Ouyang warp factor is of 0{gsM'^ /N). Terms of 0{M^) in ( [2.45 ) come from products of gsN'P /N 
with F and since F ~ 0{M) < 1, the 0{M^) < ©(M^) can be ignored. The above arguments 
could be made more precise by introducing in the estimates, but the basic idea remains same as 
above. 

^'*Note that the third term in ( 2.45 ), because of the 9i derivative, is suppressed as g^M'^Nf. Using 



( 2.28 ) and footnote]^ this would go to zero as e^/^. Also comparing this term with g^M, the fall-of 
is g^MNf which from ( 2.29| ) goes to zero as e. Therefore from all criteria in (2.45), ~ 0{M) 
seems consistent. 
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The key point in the above argument came from Li^joi'^ ^ 1 appearing in the 
Ouyang solution, which is on a regular cone while we have a deformed cone. How 
can we use the form of as given by (|2.46|) for the case of a deformed cone? The 



answer lies in the fact that for large radial distances, the deformed cone coincides 
with the regular cone. The Klebanov-Strassler solution in the large r regime behaves 
as the Klebanov-Tseytlin solution, i.e the warp factor for KS model becomes 



/2 



gtM'logb + gtM^og 
1 H 7^ log 



N 



(2.47) 



where h is some scale and = QgNa'"^ with = QsMHog b,rQ = hr^,. The above 
expansion shows that the KS warp factor in the deformed cone can really coincide 
with the Klebanov-Tseytlin solution. Once back- reactions of the flavor D7 branes are 
taken into account, KS solution in the deformed cone background will take the form 
of the Ouyang solution. We can of course choose log 6^1 such that L'^/a'^ ^ 1, 
so our argument that M/N <^ 1 holds even if we started with KS solution and not 
the Ouyang solution^^. Hence it is justified to use the Ouyang solution even for the 
deformed cone. 

Also note that, although there were no D3 branes in the KS solution, an effective 
= QsMHog b reappears in the warp factor of KS model in the large r region. This 
N can be identified with the appearing in the Ouyang solution which also justifies 
using the Ouyang solution on the deformed cone background for large r region. For 
small radial distances, we cannot use the as given in ( |2.46| ) — hence the non- 
extremal solutions we consider are only valid for large radial distances. This also 
means, we are considering large horizon Vh and the geometry is dual to the high 
temperature regime of the gauge theory. A conclusion that is consistent with our 
earlier works. 

2.1.2 Analysis of the full background with backreactions 

Once the behavior of F and the suppression orders for various terms are laid out, we 
are ready to tackle the backreactions to order QsNf and g.sM'^ /N. We start from the 
equation of motion for r given in the following way: 



(2.48) 



However, the underlying F-theory picture [|T7[ on which we based our solution 



dictates that dr ~ 0{gsNf) and therefore we will ignore terms of 0{g^^Nj). So the 



^Incidentally, using we would require b to go to infinity as exp (e 



-9/2\ 
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precise from of r will not appear in any of the equations ( p.45|) , ( p.l5 ), ( 2.33 ) and 

(PD. 



Thus with our ansatz for the metric (p.l| ), ( |2.3| ) and choice of fluxes, we have four 
equations (|2.45|) , ( p.l5| ), ( p.33|) and ( p.35p that we need to solve and three unknown 



functions h^,G and F. However, it is more convenient to write ~ A^L'^/r'^ and 
then from (|2.45|) one readily sees that 



~ 0{M/N) < 1, with F ~ 0{M/N) + 0{gsM^/N) + 0{g^^M^Nf/N) < 1 

(2.49) 

and so the third term in F is even more suppressed. Now what can we say about G? 
As already pointed out, G ~ 0{M, ggNf). But using the form of F as given above 
in ( p.49| ), one readily gets from expanding ( p.l5|) , that 



G = 0{F) ~ 0{M/N) + 0{gM^/N) + 0{glM^Nf/N) < 1 (2.50) 

Thus it is reasonable to consider only up to linear order terms in A^, G and F. But 
( ^.351 ) is a trivial equation up to linear order (see Appendix ^) and hence the only 
non-trivial equations we are solving are ( p.45|) , ( p.l5|) and ( p.33|) . Thus we have a 
system of three equations and three functions A^^G and F — which can be easily 
solved. 

Note that once the above three equations are solved, the corrections from all the 
other Einstein equations are automatically suppressed, as long as we are in the range 
( p.31| ), and the precise functional form of the axion-dilaton field r and the non-ISD 



three- form flux G3 do not influence the four equations upto linear order in A^^G 
and F. This is because ( |2.45| ) is obtained from ( |^.l(j| ) which is identical to ( |2.13| ) 
(up to linear order in A^,G and F) which in turn is obtained by tracing Einstein 
equations in the Minkowski directions. On the other hand, ( p.33|) is obtained from 
tracing the Einstein equations in the internal directions. Hence a solution to ( p.45| ) 
and ( p. 331) along with the background Ouyang warp factor and three form fluxes 



G3 minimizes the action ( |2.5|) where only Ricci scalar and the flux strength appear 
for the radial range ( |2.31| ). Thus solving ( p.45| ) and ( |2.33|) really means putting the 



action on shell which guarantees that individual Einstein equations change the metric 
only to order r^/r^ as depicted in ( p.3| ). 

The form of the solutions to the three equations along with the boundary condi- 
tions that dictate the behavior of the warp factor A, B near the horizon is discussed 
in Appendix Here we only quote the functional form of the solutions 

= ^ [Aq + Ai log r + A2 logV) 

-4 -4 

e^^ = (7 = 1 - ^ + G = 1 - ^ + (70 + log r + logV 

F = Fo + Fi log r + F2 logV (2.51) 



-17- 



where Ai, gi, Fi for z = 0, 1,2 are in general functions of r and the internal coordi- 
nates 6j, (f)j, ip, with j = 1,2. In Appendix ^ we have worked out the simplest case 
where Ai, Qi, Fi are assumed to be functions of r only by neglecting 0{gsNf) terms^^. 
This is a reasonable assumption for small number of flavors. Furthermore, the ther- 
modynamics of the field theory is dictated by the behaviour of the dual geometry 
near the black hole horizon ( [2.311 ). If we keep all the seven branes away from the 
black hole, we can ignore running of r near the black hole. On the other hand, for 
constant r we expect a Klebanov-Strassler type solution which essentially means the 
warp factors A, B and squashing factor F are only functions of r. Hence, as long as 
we are dealing with the light degrees of freedom that arise from the deformed cone 
ignoring the back reaction of seven branes, we can neglect the contributions from the 
seven branes far away from the black hole and consider the solution in ( |2.51D to be 
functions of r only. 

To account for the heavy quarks, we have to include 0{gsNf) terms but our 
ansatz ( |2.51| ) remains the same with the understanding that now Ai,gi,Fi are addi- 
tionally funtions of the internal coordinates. Interestingly, however, to analyze the 
melting of the heavy quarkonium states, we can consider string world sheets that are 
fixed in the internal directions which results in evaluating the warp factors A, B only 
for fixed values of the angles 9j, ip. This means our above analysis would suffice. 
Hence, even for the study of linear confinement and melting of heavy QQ pairs, it 
is sufficient enough to treat the solutions in ( p.51|) as being functions of the radial 



coordinate only (see [|29|, ^ for related works in this direction). 

In Figures 1, 2 and 3 we have plotted g{u), Aq{u) and -Fo(^) where u = r/fh 
using the numerical solutions to equations (|2.45|) , (|2.15|) and ( |2.33| ). As discussed in 



Appendix |C|, at the lowest order of perturbation, only keeping upto linear order 
terms in ggM'^/N, equations ( |2.45D , ( |2.15| ) and ( |2.33D drastically simplify. We obtain 



a solution with only Ao,go and Fq non trivial while Ai = A2 = gi = g2 = Fi = 
F2 = 0. For the plots, we have chosen gsM'^/N = n/S and the following boundary 
conditions^'' 

Ao{oo) = 0, Aq{oo) = 
5(0(00) = 0, 5-0(00) = 

Fo(oo) = 0, F(;(oo) = (2.52) 

Note that 5(1. 04) ~ 0, indicating that the horizon has shifted from the AdS black 
hole value of and we have obtained a larger black hole with horizon Vh ~ 1.04 fh. 
Our numerical results shows that go,Ao, \Fo\ < 1 which validates our perturbative 



i^It should also be clear that A, ~ 0{M/N) from ( |2.49D . 



^^Let us assume, for simplicity and for performing the numerical analysis, QsM"^ /N = 7r/3 to be 
the smallest scale in the theory (instead of M/N that we took earlier). Then the argument used 



earlier in (2.31) will imply that we should trust our result for r > 0.98r/j. 
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analysis. The fact that the black hole is of larger size than the AdS limit is consistent 
with the underlying gauge theory structure^*. The presence of the fractional branes 
has increased the effective mass of the black hole. In fact, the black hole entropy is 
larger than the corresponding AdS limit since Ao(r/j) > and using Walds formula, 
one readily gets that s/T^ ~ N^^ > N"^ where we have defined gsNes = fh^ijh)- 



8 

1.0 




Figure 1: The blackhole factor 5 as a function of u = r/r^. We have plotted g along the 
y-axis and u along the x-axis. Using above considerations, one may trust the result for 
u > 1.02. 



Finally, note that the identification of F with in ( |2.4| ) implies that the reso- 
lution parameter is given by 



+ 7r (^-^/) ""h Pi4— ^ + log sin-sm- (2.53) 



a 



47r \ A / \ ^ 

where we show the bare resolution parameter^^ in F and a?. The coefficients Pij 
are constant numbers that could be determined from ( p.51| ) and Appendix 0. The 
above representation of the resolution parameter is perfectly consistent with our 
conjecture in [^]: the resolution parameter will pick up dependence on the horizon 
radius r^. Interestingly we now have managed to get the leading order log r 
corrections to the result. 



^^Also note that the result is consistent with the first law of black hole thermodynamics which 
states that the increase in horizon radius is related to the increase in the mass of the black hole. 
The addition of five-branes have increased the effective mass of the black hole compared to the AdS 
limit. 

^^In the limit where the bare resolution parameter vanishes, which is the Klebanov-Tseytlin 
solution, we see that the gsM'^/N corrections actually make the small r regions non-singular creating 
an apparent resolution parameter proportional to the horizon radius. 
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Ao 




Figure 2: Non-extremal contribution to the warp factor given by Aq plotted as a function 
of ti = r/f/j. j4o is plotted along y-axis, and u is still along the x-axis and our analysis is 
valid for u > 1.02. 




Figure 3: The squashing factor given by the resolution function Fq as a function of 
u = r/fh. Note that Fq is always negative for all distances outside the black hole horizon. 
We plot Fq along the y-axis, and u is along x-axis, as before with validity from u > 1.02. 



However there is one issue that might be confusing the reader. From Figure 3 
we see that Fq is always negative for all values of r in the range < r < oo. Our 
identification of F with would then imply a to be a purely imaginary number. 
However surprisingly this does not create a problem. As we will show in the next 
subsection, all the fluxes etc. are completely expressed in terms of a^, so that a does 
not appear anywhere. Even terms with logarithms, for example ( 2.80 ), appear as 
log |a^|, so that < do not create any inconsistencies. This is of course shouldn't 
come as a surprise because the resolution parameter appear in the metric ( |2.3|) as 
1 + Fo and since |Fo| < 1 it shouldn't lead to any inconsistencies no matter how we 
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relate Fq to a. 

In our opinion the result that we presented above is probably the first time where 
the backreaction effects from black hole, including the resolution factor, are taken 
into account in a self-consistent way to lowest orders in QsNf and QsM'^ /N. To this 
order, as we showed above, the backreactions from fluxes and branes could be consis- 
tently ignored in the near horizon limit ( |2.31| ). One may now take this background 



and compute the IR effects for large thermal QCD. However before we go about 
studying these effects we would like to dwell, just for the sake of completeness, on the 
corrections to the Klebanov-Strassler three-form fluxes that arise from the backreac- 
tions of the black-hole, local brane sources, and the resolution parameter. Readers 
wishing to know our results may however skip the next sub-section altogether and 
proceed on with the calculations of the RG flows and the effects of the chemical 
potential. 

2.2 The three-form fluxes revisited 

In the above subsection we managed to provide a detailed derivation of the non- 
extremal limit of the Klebanov-Strassler type solution with a background warped 
resolved conifold. The ansatze that we used to solve for the fluxes was ( p.39| ) where 
we divided the three-form fluxes into two pieces: one coming from the known Ouyang 
fluxes, and the other coming from the various backreactions. The second piece, for 
both RR and NS three-form fluxes, received contributions from the bare resolution 
parameter Og and the QgM'^/N terms in addition to the 0{rh) terms. In the follow- 
ing we will not only justify this but also provide the form of the three-form fluxes 
including the above-mentioned corrections in the limit where the second squashing 
factor Q in (|2.3|) is negligible. Our analysis will also not be affected by the constraint 
( p. 31 ) that we had to impose to solve EOMs in the above subsection. In particular 



this means that the radial coordinate may take all values above Vh- 

Using the metric ( p.l| ) and ( |2.3| ) with the condition ( |2.4| ), the non-ISD RR three- 
form flux F3 = F3 — CqH^ takes the following form^°: 

F3 = (a. - ^;^) Y: (sin 9, de, Ad<P,-J2^ sin 9, d9, A 
A — — > — dr Ae^ A cot — sm 9^ dcho — > cot — sm 9i 

a \ a 

3gsM{r)Nfea ■ n ■ n ( J2 ^sr^ K . 
- 2^ 8nrn^^ '''' '''' ^ + 7^ 

a \ a 



d92\ A d(j)i A d(j)2 (2.54) 



^'^For the derivations of the three-form fluxes, the readers niay want to look up our earher papers 
|, |, |l6l where ah the necessary details are given. For example, the ISD fluxes on the resolved 
conifold are derived in [||, and their extension to the non-ISD cases are argued in ||, |l6|. In the 
following we will elaborate more on the derivations of |5| |l^ and show the consistency of the results 
presented therein. 
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with tto 



1 + ^ and is defined in the intermediate region rn 



< r < r„ 



The 



additional contributions to (|2.54|) are all proportional to powers of Vh, as they vanish 
in the ISD case. Finally, the quantity e^a) is defined in the following way: 



-(«) 



a 



+ E 



''an 



(2.55) 



with hem are functions of gsNf, M and the horizon radius r^. In a similar fashion 
P{a), S{a) etc are also defined. The other coefficients, for example Cq, ...ha would 
again be functions of QsNf and r/j, but also of the resolution factor a (including the 
internal angular coordinates). The resolution factor appears from the gravity dual 
that we considered in i.e a resolved warped deformed confold with the resolution 
factor can be viewed as a function dependent on the horizon radius Vh. We will argue 
this in details below. The coefficients ban can be represented in terms of the following 
matrix: 



b{an} = 



/ ^00 ^01 ^02 ^03 

bio bu bi2 bi3 

^20 ^21 &22 &23 

^30 ^31 ^32 ^33 

\bmO b 



(2.56) 



The elements of the matrix ban can be determined in terms of the cq, ci, C2, .. coeffi- 
cients that appears in the expansion ■;^t^. This is one reason of writing the various 
powers of r using different symbols. For example cT(q,) will have a similar expansion 
as ( p.55|) but with a different matrix. The various elements of the matrix will now 
be determined in terms of do, di, d2, .. etc as one would expect. For the first case, we 
have managed to determine Cq, upto few terms. They are represented as: 



Co 



l + 0{rh), ci = 0{rh), C2 



277(2 



logC +0{rh) (2.57) 



where a is the resolution factor and C is a parameter whose importance will become 
apparent soon. Once we know Cq,, it is not too difficult to get the relations between the 
various components of the matrix ( p.56| ). One may now show that the components 
ban satisfy: 



^00 = boi = bio = 0{rh) 

Cobo2 + Cifeii + C262O = 



2coboobo\ + C16 



10 



0(r,) 



C0601 + 2co6oo^o2 + 2ci6io6ii + C2620 = 0{rh) 



(2.5J 



- 22 - 



Following the above set of relations one may show that: 



E 



[2 -Slog {Cr)] + 0{rh,glNf 



An 



(Cr) 



(2.59) 



which is consistent with what we discussed in namely, the resolution parameter a 
can be thought of as a function of (r/^, M, gsNf), including the radial and the angular 
directions, i.e 



a\rh,M,g,Nf) = a^ + J^ 



a=l 



b^gf[M{r)rh 



la+1 



(2.60) 



with ha being functions of the angular directions so that this is consistent^^ with 
( 2.53 ). These can be determined by comparing ( p. 601 ) with ( p. 53 ) derived in the 
previous subsection (note that ( |2.53|) implies 6o = 0). The other similar factors 
appearing in the flux ( p.63| ) are given in terms of the following series expansions 
similar to (12.591) above: 



a 



1 + 



l%a\rh,M,gsNf)\og ((r) 
lSa'{rh,M,gsNf)\og ((r) 



+ 0{rh,M,g,Nf) 



+ 0{rf„M,g,N 



(2.61) 



Note also that there are squashing factors given by -fe- and -jtt- These squash- 
ing factors distort the spheres and therefore affect the fluxes on them. Its easy to 
show that these factors are given by: 



729 glN]a\rh,M,g,Nj) 



log (Cr) [2 - 31og (Cr)] + 0(r,, M, g^Nf 



327r2 



(Cr) 



81 ^,iV^a^(r/„M,^,iV^)log (Cr) 
Stt ■ (Cr)2 



(2.62) 



9,]V,a=(rk,A/,9.iV,)log ((r) 



E 



9a 

/pP{ct) 



1 + 



2 47rr2 + gg^Nja^irh, M, gsNf)[2 - 3 log (Cr)] 
36a2(r;„ M, gsNf)\og (Cr) 6A8a\r,,, M, gsNf)\og\C7 



+ 0{r,,M,giNj) 



+ 0{r^,M,glN}) 



a 



(Cr) 



(Cr) 



^""^Its not too difficult to argue for ( 2.60 ) using ( 2.45 ) and ( 2.46 ) if we say that goes as 
i^p{M/N) + OigsW^/N)] from (^.49|). Tfie first term in F is of 0{M/N) from ( ^loj) which 
means Oq is of 0{M/N). This will again be shown later in this section using a slightly different 
argument. Once this is established, comparing both sides of ( p.45| ) then easily implies ( p. 60 ). To 
lowest order then a? will be a function of (rf^, M/N, ggM"^ /N, ggNf) as shown in ( [2.53 ). 
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1 + 



^Qa\rh,M,gsNf) log ((r) 
(Cr)3 + 18a2(r,, M, (?,iV^)Cr log 



+ 0{ru,M,glNf 



^ _^ ^ ^ _^ 36a^(r;„M,^,iV;)log ((r) ^ 648a^(r,, M, ^,iV^)log^(Cr) ^ ^^^^ ^ ^^^2 



1 + 



36a^(rfe,M,(7,iV/) log ((r) 
(Cr)2 + 18a2(r,,M,^,iVy) log ((r) 



{Cry 

+ Oir,,M,g'^Nj) 



The far IR physics is then determined from ( p.63| ) and the squashing factors ( p.62| ) 
by making the replacement {(r) r to the radial coordinate. Note also that all 
the flux components are expressed in terms of and therefore the sign of can 
be directly inserted here. Considering all the above, this then gives us exactly the 
result that we had in namely: 



Fs = 2MAi (l + %^ log 



SgsMNf dr 



An 
SgsMN 



r ] A - (sin 9i dOi A — Bi sin 62 d92 A d(f)2) 
' 92 d(f)2 — B2 cot — sin 9i dc. 



A2 — A A I cot — sin 9: 



f 



Stt 



9' 



9i 



A3 sin 9i sin 92 cot — d9i + B3 cot — c?6'2 A d(f)i A d(j)2 (2.63) 



where we have taken M(r) — )■ M in the far IR, and the various coefficients Aj, Bj 
are related to ( p.59|) and (|2.62|) as: 



E 

a 



fa 



r 

Qa 



Ai, > = A2, > 



B2, 

^' ZL^ ^<5(q) 



A3 
B3 



(2.64) 



As we mentioned before, the additional contributions to ( |2.64| ) are all proportional 
to powers of Vh, as they vanish in the ISD case. Needless to say, the functional form 
for F3 is consistent with ( [2.391 ). 

The NS three-form flux H3 is now more interesting. Unlike F3, it has to be 
closed. When the resolution parameter a and M are just constants, it is easy to 
construct a closed H3. In the presence of non-constant a and M(r), finding a closed 
H3 is more non-trivial. For our case H3 is given by^^: 



^^3 = E 



6g,M(r)k„ 



1 + 



1 (cosec Y cosec 
2^ 



I) 



9gsN-f 

2-Kr 2 



A^ I sin 9i d9i A c/0i - ^ sin d92 A rf02 j + ^ 

\ a /a 



dr + J2 0{rh)da. 

i 

^glM{r)Nfl^ ( dr 



^ Ae^--de^ 



^^We correct a minor typo in pf 
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A l^cot I - cot I d9^^ + 9s^^ (&i(r) cot | d9, + ^^(r) cot | rf^^) 

3gsdM(r) / 1 Oa^^.A^f 



1 log 

















47r (ir y ' ^ r^gsNf ^2 ^ 

sin 9i dOi A d0i A dr - — 2 + 9^,A^f - — ^^^^ 

^ ^ 127r dr \ ' j ^wg.Nj ^2gsNf 

log I sin — sin — ) sin 62 d92 A (i02 A dr — — log (sin — sin — ) sin 62 d92 A d(j)2 A rfr 

\ 2 2 / 1271 dr V 2 2 / 

(2.65) 

with {ka,...,qa) being constants, cXi = (61^,0^) and 6„ = ^Jf^ where = 
{gsNf,M,rh) and {9sNf,M, Th). 

Note that in addition to the simplest 0{rh) term that we mentioned above, there 
would be more terms of the same order that vanish in the ISD limit. We might also 
wonder about terms of the form da/dr and da/dcji. From the form of ( |2.60| ) we 
see that these terms themselves are of 0{M) so to this order they could either be 
absorbed in dM{r)/dr terms or in the 0{rh) terms. The various squashing factor 
etc are now given by: 



a 

^ Pa _ ^ , ?>9sO?{rh,M, gsNf) 9gsa'^{rh,M, g^Nf) ^ ^ 



,2/^. ]\/r „ AT A 1^^ /'/•^^ Oc;no^4/^. T\/r „ at \ i^„2 



^0{r^,M,glN)) (2.66) 



12a\rh,M,g,N^) log ((r) , ,^^2^ 



Cr + 36a2(r;„M,(/,iV;) log ((r) 

To study the far IR physics, we again consider M(r) — )• M, with the various expan- 
sions in ( p.66| ) are related to A4, A5, B4 and B5 respectively. Note again that the 
resolution parameter in all the coefficients appear as a} . This then reproduces again 
the far IR result of |T6[ as well as the expected ansatze ( p.39|) , namely: 



= 'og^K^^M I 1 + log r + log sm— sm— | — 
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A- sm d6i A d(j)i - B4 sm 62 d62 A d(j)2 + — --^5 — A e^, de^ 

2\ j m \ r 2 

A ^cot ^ - Bg cot y rf^i j (2.67) 

with the necessary 0{rh) terms that vanish when the horizon radius vanishes. 

In the far IR the closure of H3 is again non-trivial because the resolution pa- 
rameter a is no longer a constant now, although M(r) — )■ M. All the informations 
of non-constant a are captured in the coefficients A4 5 and the squashing factors 
In the following section we will determine the resolution parameter to 0{M'^) 
in the IR. This means up to this order, the closure of i^s implies the following three 
conditions: 

(z) aacot ^ — -i + a2 Agcot ^ 1^ + ^(^'^' 9sNf) = (2.68) 

Z Oui Z Oui 

.... ■ n dA, 62 dA, 02 OA, 

[ii) ctism Oi — — + a2Cos 61 cot — — assm 61 cot — ^— 

0O2 2 afc*! 2 or 

9 (9B 

+ asAgcos ^1 cot ^ — ^ + ^(r;,, M^ (^.A^^) = 
/ 002 

/XX. n 01 dA4 , . ^ (9B4 , ^ ^1 9Br 

(zzz) ai B4 sm 62 cot — — «i A4 sm 6*2 ^-^ ^2 Ag cos 6*2 cot — — — 

Z oUi Oui Z OU2 

^1 9Bg ^2 9A. 

- as Ag sm 62 cot — — h ^2 cos O2 cot — — — 

2 or 2 oOi 

- as Bg sin 62 cot — ^ + 0(r^, M^ f^.A^^) = 
2 ar 

where ai, 02 and 03 are defined as: 

3gsM ( 9gsNf g,Nf . 6, . 62 

«! = 1 H log r H log sm — sm — 

2 \ 47r 2tc 2 2 

a, = MMIk, „3 . JjiMl . (2.69) 

Svrr ' 167r 2 ^ ^ 

The RR three-form fiux F3 = F3 — Co-f^s is not closed, but it satisfies the condition 
dF^ = —dCo A i^s, which is equivalent to the statement that F3 is closed. Of course 



^1 



as described in [^, the closure of F3 is only in Region 1. In Region 2 there are anti 
five-brane sources that make F3 non-closed. The closure of F3 in Region 1 implies 
the following nine conditions on the various coefficients of the three-form fluxes: 

r ^ n^^^ ■ n n dA2 a X ^^-^ ■ n 

[iv] Pi^T— sm Oi - P2 -iT^ cot — sm Oi - P2 A2 -7^ cot — sm Oi 
Or 061 2 OOi 2 

6 9 
+kai A4 sin + ka2 Ag Bg cot y cos + ka2 Ag Bg cot y + 0{rh, M^, QsNf) = 
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(9A (9B 6 

(f ) - /3i Bi — — sin 02 + /3i Ai — — sin 62 + /32 cot — sin 62 
or or 002 2 

—kai A4 B4 sin 62 ~ ka2 A5 cot — cos 62 + A;tt2 A5 cot — + 0{rh, M , QsNf) = 

9B dA. 9 OA. 9 

ivi) j3i — — cos 92 sin 9i - 132 -777— cot — sin 9i cos 9i + /32 B2 ^r;— cot — sin 9i cos 92 
or o9i 2 c/6'1 2 

— /32 A2 ^rr^ cot — sin 9i cos 6'2 + ^rr- cot — — kai A4 sin 6^1 
09i 2 c/r 2 

- ka2 A5 B5 cot y cos 6*1 + fcas A5 B5 cot y cos 612 + 0{rh, M^, s-sA^/) = 

dA. 9B 9A 9 

ivii) — f3i Bi — — sin 92 cos 9i — Bi Ai — — sin 6^2 cos 6*1 — B2 „ „ cot — sin 6^2 cos 6*1 
ar or o92 2 

-P2 B2 ^^7— cot — sm 9i cos fc'2 - P2 A2 ^^7— cot — sm 9i cos 6'2 + B3 — — cot — 
Ou2 / Cc72 or Z 

OBs 9i 6*2 
+ A3 — - — cot — + ka2 A5 cos 9i cot — — kai A4 B4 sin ^^2 
Or 2 2 

- A;a2 A5 cot -j cos ^2 + 0{rh, M\ g^Nj) = 

/3i — — cos 92 + /^i Ai — — sm e'a cos e*! - fi^ B3 — — cot — + -ttj— cot — 
C'C72 ool Ooi Z O02 Z 

(9B 9 9 9 

-Ps A3 ^ cot 4 - ^"3 As cot ^ sin + kas A5 Bg cot sin ^2 + C(r;„ M\ QsNf) = 
O't'i z z z 

dA 9 

(ix) (3i sin 9i + ka^ A5 cot ^ sin 9^ + 0(r;„ M^, ^,iV^) = 

OU2 z 

(9B 9 

(x) /3i Ai — i sin 92 + A;a3 A5 Bg cot sin ^2 + 0{rh, M\ g^Nj) = 

Oui I 

dA 9 9 9 

(xi) - 132 -^T- cot — sin ^2 - ^"2 Ag Bg cot — cos 6*2 + ka2 Ag Bg cot — 

o9\ Z 2, 2, 

+ O{r^,,M\g.,Nf) = 

(xii) - (32 A2 cot sin 6^1 + (32 cot ^ sin 6^2 - /ca2 Ag cot 

Ou2 Z O02 Z Z 

Q 

- ka2 Ag cot ^ COS 9i + 0{rh, M\ g^Nf) = (2.70) 
where we have already defined ak, An and Bm- The (3i are now defined as: 

ft.M(l + ^log.), (2.71) 

Note that in (|2.68|) and (|2.7CI|) we have separated the 0{rh, M^, ggNf) corrections from 
the resolution parameter in the fiuxes. We may also absorb these corrections to the 
resolution parameter and write the three-form fiuxes completely in terms of 0{M^) 
corrections to the 0{M) terms in the original Ouyang solution. This may also be 
interpreted as though every fiux components sees a c?Zj^erent resolution parameter a|. 
Note also that we don't have an (9(M^) corrections to the Ouyang three-form fiuxes. 
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However from (|2.60|) we do expect an 0{M) term for a^, unless of course 60 = or 
ao is proportional to M. Comparing (|2.6CI|) with (|2.45|) , ( p.46|) and (p.53|) may imply 
bo = 0. Additionally, the scenario with Oq being of 0{M/N), is more likely, as in the 
absence of wrapped D5-branes the gauge theory is conformal with the gravity dual 
given by AdS^ x T^'^ [^]. Only in the presence of wrapped D5-branes the gravity 
dual becomes a resolved warped-deformed conifold so the resolution parameter qq 
should depend on M. This will be consistent with ( |2.45| ) and ( |2.49| ) as discussed 



earlier. In either case, it is clear that the fluxes that we take contribute the 0{M^) 
terms to the original Ouyang fluxes. This helped us to get a consistent background 
in the presence of fluxes and a black-hole as we saw in the previous subsection. A 
more elaborate study will be delegated to |T^ . 

Before we end this section, let us also see how the squashing factors in the three- 
form fluxes behave in the light of the result ( p.53| ). To 0{gsM'^ /N) the resolution 
parameter is only a function of the radial coordinate r. This means that An, Bm 
can be written as functions of r to this order satisfying the closure conditions ( [^.6(^ ) 
and ( p^.TOI ). For example, combining (^3^) and ( |2.68| ), A5 takes the following integral 
form: 

A5 = i ^ / dr^ (2.72) 

sin 61 cot Y J '^3 

where is a function of the angular 9^ and the radial r variable such that = 0. 
Similarly: 

sin 9i cot % f dr c^a'^^ir) 

B5 = f T , , (2.73) 

sin 62 cot Y J (^5<^3 v) 

where again C5 is like 05 discussed above. Once (05, C5) are determined the two integral 
forms (|2.72|) and ( p.73D not only satisfy the closure conditions ( p.68| ) but also the 



necessary EOM. These two integral forms are also consistent with the conditions 
(ix) and (x) of ( [^.7U| ) because 03 = — ^^\q^ ^ is a constant. Note however that, 
to this order, the integral form for (A4, B4) cannot be determined by this method, 
although we will know (A4, B4) in terms of the resolution parameter upto the 
0{rh, , gsNf) corrections. 

On the other hand both (Ai, Bi) do have an integral representation if we assume 
that (A4, B4) have some integral representation (which in turn will be determined 
in a different way from the one that we have followed here). If this is the case, then: 

dr f ai , A , ; A -n ^ ^1 ^ /I , ; * n 2 ^1 



Ai = / — - — — + kai A4 + ka2 A5 B5 cot — cot 9i + ka2 A5 B5 esc — 

Pi \sm Oi 2 2 

^ ^1 . r . ^ . u . . O2 „ , . o U2 



Bi = Ai / o :r + kai A4 B4 + ka2 A5 cot — cot 62 — ka2 A5 esc — 

l3iAi \sm 02 2 2 



(2.74) 
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where (ai, bi) are functions of r and 6i such that = = in the same sense as 
mentioned earher for the other cases. The other squashing factors in ( p.70|) do not 
however have such simpler integral forms. 

Another interesting thing to note is that from condition (xii) of ( ^.701 ) we might 
get a simpler form for A5, namely: 

A5 = ^ (2.75) 

2A;a3Cot %(1 + COS ^1) ^ ^ 



This may seem to be different from ( p. 72 ) that we derived earlier. This is however 



not the case because the coefficient 07 is related to 05 in the following way: 

«5 = 7:U2r + ^7J2nr: (2-76) 



k \d9i dOidr 

One may also cook up somewhat similar relation for 65 and 67 from condition {xi) of 
( [^.7UD as above. The final result will again be consistent with what we got earlier, 
establishing the fact that the system is well defined with the given set of boundary 
conditions. Therefore these analyses complete the side of story that we expected 



from H] and |16| in a satisfactory manner. 



2.3 The behavior of the coupling constants 

In the above sub-section we computed the background more or less exactly upto 
order . To this order we see that the corrections to the resolution parameter 
is only functions of r, the radial coordinate. If we go beyond this order, the angular 
dependences start showing up. 

One other thing along the same line would be to study the behavior of the two 
gauge coupling constants of the boundary theory. For example a crucial question 
would be to ask how the RG fiows of the coupling constants change when the thermal 
effects are turned on. In the literature there have been many confusing statements 
on this. In the following we will argue that the RG flows or more appropriately the 
thermal beta functions should be properly interpreted and the correct picture, in our 
opinion, is that the thermal beta functions do not change, but the coupling constants 
themselves get renormalised. Let us elaborate the story below. 

Once we know the NS 5-field B2 and the string coupling e* then it is easy to 
determine the gauge couplings at the UV of the dual gauge theory. The resulting 
relations are: 




(2.77) 
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Now note that when is a constant the string couphng and the B2 field were 
obtained in ISll as 



s = log(r +9a r ) log sm — sm — 

gs Sn 2n \ 2 2 

B2 = (bi{r) cot ^dOi + 62 (r) cot ^(i^s) A 



+ 



^ ^(1 + log(r2 + 9a2)) log ( sin ^ sin ^ ) + h^{r) 
An V 2 2 / 



sin6'i(i6'i A 



^^^^(-36a2 + 9r^ + 16r^ logr + logfr^ + Qa^)) log ( sin ^ sin ^ ) + 64 (r) 
. 127rr^ V 2 2 



X sin 62(162 A (i</)2 



(2.78) 



where we have shown the dependences in i?2- The other dependences come 
from the implicit coefficients bi{r). These dependences in the first two coefficients 
bi{r),b2{r) are given by^^: 



^i(^) = ^^^^f (18a' + (16r2 - 720^) logr + (r^ + 9a^) \og{r^ + 9a^)) 



(2.79) 



Similarly, other dependences on the resolution parameter would come from the (63, 64) 
coefficients. We have determined these coefficients in terms of first-order differential 
equations. They are now given by: 



glMN 



f 



Kir) 



SgsMir"^ + 6a2) g'^.MN^ 



,2 I oc„2 



36a' - 18a' log I a' I +34r' logr (2.80) 



;iOr' + 81a') log(r' + 9a' 



18a'-18(r' + 6a')log|a' 



+ (34r' + 36a') logr + (lOr' + 63a') log(r' + 9a' 



As all the above coefficients are given in terms of the resolution parameter derived in 
( p.53|) , they should then be functions of r/j and other radial and angular variables'^. 
This means that the exact coupling should be determined in terms of the NS 5-field 
B2 that is of the form, upto the order that we had in ( |2.53| ): 



B2 = B2{r,6,) + {gsM)-{gsNf) 



9sM^ 
N 



B2{r) + {gsM) ■ {g,Nf) 



9sM^ 
N 



C2(r,^,) 



^■^The bi(r) here are the same bi{r) that we encountered in ( 2.65 ) in the Umit where the fluxes 
become ISD. 

^^Note the appearance of log |a^| so that < will not lead to any inconsistencies, as we 



explained earlier. 
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(2.81) 



and not B2 that we mentioned in ( 2.78 ). Also we expect dB2 = with given as 
in ( p. 651) . However, notice that the extra terms in (|2.81|) are suppressed by gsMNf 
over and above the ^"j^f suppression. Therefore if we follow the limit shown in (|2.28|) , 
we can easily infer that this makes ( |2.81 ) and the resolution parameter (|2.53|) to 
have the following e expansion 



B2 



,9/2 



ri + e 



9/2 



(2.82) 



that basically tells us that B2 ~ B2 which is consistent with our assumptions in [^]. 
To this limit then the running of the couplings do not change, as one would have 
expected. This gives us: 

3iV; 



d 



9 log A 

d 

Slog A 



47r2 47r2 



91 
47r2 

11 



9l 
47r2 



4 



9l 



3M 



r2 + 9a2 

3^ 
47r 



log(r2 + 9a^) 



(2.83) 



where A is the energy scale in the gauge theory side. In fact the LHS of ( |2.83D is 
written in terms of gauge theory variables whereas the RHS is written in terms of 
gravity variables. 

One important question now is to ask what happens when we consider 
corrections. Clearly now we need to consider the corrections to ( p.53| ). What does 
this imply for the running couplings? Saying that the coupling constants run at a 
different rate would probably not be a meaningful statement. The correct thing to 
say at this stage would be to allow for new effective, couplings 'gx and 5^2 that again flow 
at the same rate as before, i.e the effective couplings have the same beta functions 
as the original theory. Alternatively this means that between the two scales, energy 
and temperature, we fix the energy scale and define effective couplings for any given 
temperature. Once we change the energy scale, these couplings should run exactly 
as before. 

A way to see this would be the following illustrative example. Imagine the 
complete corrections to (|2.83|) , from the changes in the fluxes and the resolution 
parameters, may be represented as: 




3A^/ /r^ + 6a^ 



"2 + 9a2 



(2.84) 



n=l 



(9 log A 



3Af 



47r 



log(r2 + 9a^) + 



m=l 



where J-'n(r, 6i) and Qn{f, di) capture all the corrections. For simplicity we have con- 
sidered the corrections to depend only on (r, Oi). Of course more generic corrections 
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could also be entertained here but it would only make the analysis involved with- 
out changing the underlying physics. Therefore here we will stick to the simplest 
scenario. 

The above corrections to our earlier set of equations can now be re-arranged 
to redefine two new set of couplings gi and g2 that are related to gi and g2 in the 
following way: 




(2.85) 



where we have taken an average over the angular directions so that the couplings 
are defined only in terms of r/j and r (or temperature and the energy scale, in the 
language of gauge theory). With this, one can now easily see that the two effective 
couplings gi and g2 flow exactly as ( p.83| ) and therefore the theory has the same 
behavior, which in turn means the same renormalisation group flows, in terms of 
these couplings. This would probably be the right way to analyse thermal beta 
functions. 



2.4 Short detour on dualities and dipole deformations 

Our final aim of this section would be to take a short detour and study the effect 
of the dipole deformations on the flavor seven-branes in the gravity picture. This 
dipole deformation, since it affects the seven-branes, should also have some effect on 
the fundamental quarks in the gauge theory. We will make some speculations how 
the dipole deformations effect the far IR picture. 

Our starting assumption would be that the solutions presented in the earlier 
subsections have isometries along 4'i,4>2 and directions. This in particular means 
that the coefficients appearing in ( p.51| ) i.e {Ai, Fi,Gi, gi) are all functions of {r,9i) 
only and not of {(pi^ip). This is not a strong assumption as we saw earlier that even 
to 0{ggM^Nf/N) the {4>i,ip) dependences do not show up. It could be that the 
background retains its isometry along (0j, ip) directions to all orders in ggNf and 
gsM'^/N, but we haven't shown this here. 

Before moving ahead let us clarify a point here. Dipole (or non-commutative) 
deformations can be studied in two possible ways. In the conformal case, one takes 
the D3-brane metric written in terms of its harmonic functions, and then use TsT (T- 
duality, followed by a shift s, and then another T-duality) to generate new solution. 
The new solution is still given in terms of D3-branes and harmonic functions, but now 
there is a background -Bns field. One then takes the near horizon limit to determine 
the gravity dual of this scenario. The gravity dual has no D3-branes, but both as 
well as = dB^sss fluxes are still present. The near horizon geometry do not change 
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the internal metric too much, and therefore analysis on both sides of the story is 
somewhat similar. 

The above criteria changes quite a bit once we go to the non-conformal case. The 
gravity dual is not simply given by taking the near-horizon limits of the D3 and the 
wrapped D5-branes. To avoid naked singularities of the Klebanov-Tseytlin form, one 
now has to deform the internal space also. This means making a TsT transformation 
on the brane side, one may not necessarily get the full gravity dual picture easily. 
This is also clear in the geometric transition set-up, whose supergravity solution is 
developed in p2|, 0. So we could do TsT transformations on two sides of the picture, 
leading to two possible different interpretations. 

Thus, once we have solutions for both sides, namely the gauge-theory and the 
gravity sides, we can use TsT transformations to deform them into various different 
solutions. In this paper we will not consider the dipole (or non-commutative) defor- 
mations on the gauge-theory side of the story^^, but concentrate only on the gravity 
side. This means, given the background metric ( |2.1| ) with fluxes, five-branes and 
seven-branes, the TsT transformed backgrounds will be related to some interesting 
deformations of the four- dimensional thermal gauge theories. These deformations 
can be classified to fall we can allow four kinds of transformations. They are listed 
as foUows^^: 

• T-dualize along one space direction say X3 then shift along another space direction 
say X2 mixing {x2,X3) and then T-dualize back along X3 direction. 

• T-dualize along 0:3 and then shift^^ along one of the internal directions that are 
isometrics of the background, namely along 0i,02 or ip directions'^, and then T- 
dualize back along X3 direction. 

• T-dualize, shift and then T-dualize along internal directions. The shift will mix 
two of the internal directions in some appropriate way. 

The first operation will lead to a non-commutative gauge theory on the D7-branes 
with [0:2,3^3] = "^-823 as our algebra. The second one is more interesting. T-dualizing 
along X3 but making a shift on the directions along which the D7-branes are oriented 
i.e along (pi, 61 and ip (recall that the D7-branes wrap the two-sphere parametrised by 
(^1, 0i) and are spread along (r, ip) directions) will lead again to a non-commutative 
gauge theory on the D7-branes. On the other hand, if we make a shift along the 
orthogonal direction parametrised by 02, then the theory on the D7-branes will 

^^The dipole deformations on the gauge theory side, at least in the far IR and in the local case, 
has been discussed earlier in The readers may refer to those papers for more details on the 

multiply allowed dipole deformations. 

^^We will use {xo,xi,X2,X3) as a convenient reparametrization of {t,x,y,z) used earlier. The 
former will be more convenient for the next couple of sections, 
gain mixing X3 with one of the internal directions. 

^*For simplicity we will only consider the isometry directions. 
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be a dipole gauge theory. For the last one is T-duahzing and shifting along the 
directions of the D7-branes will again lead to non- commutative theory on the D7- 
branes, whereas if we shift along 02 but T-dualise along the D7-brane directions, we 
will get dipole theory on the D7-branes. 

To analyse these in case-by-case basis, let us study the first kind of deformation 
first. We choose the shift to be 

X2 I—)- -|- sin6'a;3, 0:3 1— )■ cos6'x3 (2.86) 

cost' 

After the series of transformations discussed above, i.e TsT, the metric (|2.1|) , be- 
comes : 

ds^ = [-gidxl + dxl + J{dxl + dxD] + \/h{g2^dr^ + dMj) (2.87) 

with the Lorentz breaking deformations along {x'^,x^) directions specified by J. 
There is also a background 5ns field that accounts for the non-commutativity. Both 
J and the i?Ns field are defined as: 

J-^ = sin^ 9h-^ + cos^ 9, B23 = tan^/i^V (2.88) 

The metric has the same form as in and the gauge theory on the D7-branes 
become non-commutative in the and x^ directions. 

For the second kind of deformation we follow similar procedure as above except 
that now we shift along ip direction and T-dualise along x^ direction. The resulting 
metric is 

ds^ = (-gidxl + dxl + dxl + 7. 2nl 2 • 2. ^^3) (2-89) 

^yh V 9 cos^ + r'' sm J 



\2 



r'^{dip + cos Oidcpi + cos 6^2^02, _^ 



9 cos^ 6' + sin 6 

where note that the x^ and the ip circle is non-trivially warped. The dotted terms are 
unchanged from the original metric (2TT). However the -Bns field now is non-trivial 
because of the ip fibration structure: 

B = T: • 2n dX3 A [dtp + COS + COS ^2#2) (2.90) 

9 cos^ + sm u 

The scenario now is interesting because we have three components of the -Bns field, 
with two of the components B^^ and B^^-^ parallel to the D7-branes and one com- 
ponent -83,^2 having one leg orthogonal to the D7-branes. Existence of these three 
components would lead to a complicated theory on the D7-branes that in some limit 



^^For this section we will ignore the 0{gsA''P /N ^r'^/r^) corrections to the internal metric (2.3). 
A more precise result will not change the physics to the order that we are studying here. 
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may be considered as a combination of both dipole and non- commutative deforma- 
tions of the world-volume theory on the D7-branes. 

As an example for the third kind of deformation^^ we first T-dualize along ip 
direction and then shift along 02 direction. The resulting metric, keeping both the 
squashing factors {F,Q) in (|2.1|), will look like 

ds^ = -^{-gidt^ + dx^ + dy^ + dz^) + \fh r^J{dilj + cos6'i cos6d(j)i + cos6'2C?02)^ 

+ ^r2j(l + F)(l + g) sinOldct^l + ...] (2.91) 

with the dotted terms being the terms unchanged from the original metric As 
expected, note that both the fibration structure as well as the 02 directions get 
warped by J. As before there is also a _Bns field. Both J and 5ns are given by: 

5^02 = ^ta.ne ?,m^d2{l + F){l + g)J 

j-i = Ocos^^ + -/ir2(l + F)(l + ^)sin2 02sin2 (2.92) 
6 

This solution shows that the gauge theory on the D7-branes has become a non-local 
dipole theory^^. 

All these solutions generated from our TsT duality operations lead to new gauge 
theories on the D7-branes. As mentioned earlier, it is not clear to us whether these 
deformations are the corresponding gravity duals of the respective deformations on 
the gauge theory side of the picture. One thing however is clear: due to the dipole de- 
formations on the D7-branes, the KK masses of the fluctuations are different from the 
original theory. In fact the dipole deformations (along appropriate directions) tend 
to make the KK states heavier |^ . Therefore we would expect operator dimensions 
on the field theory side to also change accordingly. 

The TsT duality operation doesn't change the warp factor nor the BH factor. 



Naively applying the criteria from |16| one would think this doesn't change the 
thermal behavior of the theory. However, notice that now there is an extra non- 
constant -Bns field that cannot be gauged away. This means when we write down 



■^°Note that we cannot construct another theory by shifting along 62 direction and then T-duahng 
along 02 direction because 62 is not an isometric direction. A shift along the non-isometry directions, 
like the 62 direction, will destroy the existing isometry directions making the T-duality opertions 
highly non-trivial. 

■^^On the other hand, the shifts and the duality directions that we choose are not the most generic 
ones. We can make numerous other shifts. One simple example could be as follows: we T-dualize 
along space direction ^3, then shift as z 1— > z + A0|/2 and finally T-dualize back to generate a 
non-trivial background with the metric 

ds^ = -^{-gidf + dx^ + dv'^ +dz^) + \/h{g:2^dr'^ + dMl + \^ eldel) 
and a -Bns held, -B302 = -^^2- This would generate dipole deformation on the D7-branes. 
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the Nambo-Goto action for the strings, the effect of the i?NS held can no longer 
be ignored and it'll definitely change the criteria of the confinement /deconfinement 
transition studied in |16|, This shouldn't be surprising because one would expect 
the thermal behavior of the dipole deformed quarks to be different from the un- 
deformed ones. 



3. Background chemical potential and backreactions 

After carefully working out all the details about the background, it is now time to 



study various applications. In earlier papers |T^, |20|, [T8| we managed to study 
aspects of phase transitions, quarkonium meltings, viscosity to entropy ratio and 
many other things despite lacking a precise knowledge of the background. Of course 
our ignorance about the precise details of the background, first proposed in p , wasn't 
quite a handicap as the knowledge of the existence of an UV completion was enough 
to get results for many of the above mentioned computations. For those computations 
that needed specific details, for example the question of how much we deviate from 
the well-known viscosity to entropy bound, were left in their functional forms. Once 
precise background informations become available, these functional forms may be 
replaced by their exact values. 

In the following sections we will however not re-address those questions here. In- 
stead we will ask what happens if we switch on a chemical potential in our theory. In 
fact, even before we study the consequences of switching on a chemical potential, we 
should ask how to generate a chemical potential in out set-up. This is unfortunately 
not so simple as the AdS case. Part of the reason is that, a-priori there are two 
known ways by which we could generate chemical potential: one via duality chasing, 
and the other via the D7-brane gauge fields. Both give different result, and we will 
argue below that the former, i.e the duality chasing method, is inherently flawed due 
to certain subtleties associated with the curvature of the background. The latter will 
turn out to be more useful. 

The second reason that makes the story a little more involved has to do with the 
underlying RG flows in the dual gauge theory. Due to this, we need to understand 
the effect of the chemical potential at a given renormalization scale. In the gravity 
side, this means that we have to study the background at a chosen radial coordinate. 
We will discuss this below, and also point out the subtleties with the duality chasing 
idea that make it difficult to use this as a viable method of generating chemical 
potential in our theory. But first let us give a brief discussion of the concept of 
chemical potential in thermal field theories. 

3.1 Chemical potential in thermal field theories 

Chemical potential is associated with conserved current. Once there is a conserved 
charge Q = J j^d^x, there is a chemical potential fi conjugate to Q, and the grand 
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partition function for the thermal field becomes 



y^V^exp d*x{£E-l^f) (3.1) 



where H is the Hamilton and is the Lagrangian density after wick rotation. We 
have omitted other gauge group indices for simplification. 

In the field theory the only field naturally couples to j° is so the chemical 
potential in the thermal field theory is defined as the time component of the gauge 
field Aq. Loosely speaking, according to the AdS/CFT correspondence this is the 
value of the bulk gauge field At{r) at boundary, i.e. radial infinity. If we have a gauge 
field At in the bulk its value at infinity will introduce a chemical potential associated 
with some conserved charge for the gauge theory living on the boundary. In the 
following section we will first try to generate such bulk field by duality chasing. 

A more interesting thing to do is to study the chemical potential associated 
with the quark number. In our case we have M coincident fractional D3 branes 
(in the far IR) and Nf coincident D7 branes so the SU{M) gauge theory has a 
global U{Nf) ~ SU(Nf) x U(l) symmetry. The fundamental matter ijj and 'ijj, 
which are quarks, transform under the global U{1) symmetry with charges 1 and 
— 1 respectively. Thus the U{1) charge counts the number of quark numbers Uq. In 
terms of partition function it can be written as 

Z = exp(-/3W) = exp [-I3{H - //n^)] 

where W is the Gibbs free energy which satisfies 

f = -n, (3.2) 

Prom the string point of view, the Gibbs free energy is proportional to the on shell 
D7 brane action, 

m = / d'^x dr dQsSCE = / d'^x dr dQs , drSAo(r) (3.3) 

J TO J TO odrAo{r) 

where Tq is some scale that will be made precise later and is the internal three-cycle 
along which D7 branes are extended. On the other hand 



'SdrAoir) 



The reason is the following. The U{1) field on the D7-brane worldvolume can be 
thought of as arising from fundamental strings dissolved into the D7-branes. Thus 
the density of the quarks is precisely the density of the strings. Since the fundamental 
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strings sources the NS two-form i?Ns = -Bor, so its density can be determined from 
the local charge density of the -BNS-field. Furthermore, the gauge invariance requires 
that the D7-brane action only involves the combination B^r + 27r/gFor so that 

ddrAo{r) J dBor 



Plugging ( p^) into ( pT3|) we find 

6W = -n,[6Ao{oo)-6Ao{ro)] (3.5) 

since Ao(ro) always vanishes, comparing ( p.5| ) and ( |3.2[ ) we get /i = y4o(oo). We will 
calculate this chemical potential in the next subsection after we clarify the subtleties 
associated with duality chasing method to generate chemical potential in our model. 

3.2 Chemical potential using duality chasing 

The duality chasing idea is somewhat similar to the technique that we used to gener- 
ate non- commutative and dipole theory on the D7-branes. Instead of using shift, we 
will apply boost and then T-dualize. We call this technique as TbT where b stands for 
boost. The boosting mixes Xq with the internal coordinate, say y™, to generate a -Bns 
field of the form -Bom- This means from the five-dimesional point-of-view this would 
indeed be a gauge field Aq = -Bom- Since the duality chasing preserves^^ EOMs, our 
method should give us a background with a vector potential Aq. Unfortunately, as 
we show below, this fails in a rather subtle way. 

Our starting point would be to consider the metric and the fiuxes presented in 
section 2 but ignoring the 0{gsM'^ /N,r'^/r^) corrections (this means we are in the 
range ( |2.31| )). Since our background is divided into three Regions, namely 1,2 and 
3, the easiest way would be to keep all the warp factors undetermined so that after 
TbT we can get the results for the three Regions simultaneously. This is easier said 
than done because the background also has i?NS fields that start decaying faster as 



we approach Region 3 ||T6| 



To avoid these subtleties, let us study the effect of the TbT operations only in 
Region 1 . The 5ns field in Region 1 will be of the form ( |2.81| ) . However due to the 
scaling arguments given in ( p.82| ), the 5ns field components take the form ( p.78| ). 
We can simplify this a little further by making the components independent of 02 
direction with a gauge transformation: 

Enow = B2 + dCi, with Ci= B^^i d(f)2 (3.6) 



where B2 is given in ( |2.78|) and -Bnow is the -Bns field that will enter the duality 



relations. Taking this into account, after TbT the metric will be: 

ds^ = —= ( r, "^^^ ^ — dt^ + (ixfnq I + ^^^dr^ 

^\ A cosh^ /3 - ^ sinh^ /3 (71 



32 



Both sides are weakly coupled, and so no non-perturbative effects could enter here. 
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cos^ e^A I 1 



+ 



6 



[de\ + sin^ ^1 



v^r2(l + F)(l + e;) 



9 



1- 



- cos 6'i cos^ Q2 



1 



X 



1 
3 



- (^^ sin^ ^2 + - sin^ ^2 cos^ 
^ cosh /3 



^cosh'/3- ^sinh'/3 



X 



y^r2(l + F)(l + e;) 



cos 6^2 sin^ Q2Adil)d^2 + cos 6*1 cos Q^f^ 



/ir^^cosh/3 cos^^2(l — cosh/3) 
"81 ^cosh^/?- ^sinh^/? 



{dip + cos 6* 



sin^ 6*2 c 



where more details are given in Appendix |A[ The background now indeed has a 
gauge flux Aq coming from ( |A.9D alongwith a dilaton 0. The various components are 
given by: 



cos 9o 



(3.7) 



A 



An 



y/hr^ [2cos^| + 3(l + F)(l + ^) sin 9o 



6- - In 
2 

Bn^-, d(h 



cosh^ (5 — Qi sinh^ /3 
1 + 5o^2#2 + 5o^# + Boe^dOi 



(3. 



with -Bo<}ii, -Bo</)2) and i?o6»i are the four-components given in ( |A.9| ). There are also 
the RR fields that we haven't discussed here. They do not mix with the metric and 
the i?NS fields and therefore do not change the chemical potential. This is because 
the gauge field in ( |3.8|) come from the cross-terms in the metric generated by the 
boost operation and not from the fluxes. 

Unfortunately the solution presented in ( |3.8|) and ( p.7|) has numerous issues. We 
notice that the dilaton becomes imaginary and gu becomes infinity and changes sign 
at r = ri, where 



1 



ri 



3^6 coth 13 



9i V6cos2^2 + 9(l + F)(l + ^)sin^^2 



(3.9) 



and becomes negative when r > ri. We will asume that ri >> {N/MY^^rh- When 
dilaton becomes imaginary that means the theory is not well defined in this region 
and signals a serious problem for the TbT done for this metric. Actually after the 
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boost in type IIA the metric already have gu and g4,2if>2 changing sign at r = ri 
which means at that point t becomes spacehke while 02 becomes timelike. In ^1 



the authors discussed the dimensional reduction on a timelike direction, but for our 
case its not clear whether the dimensional reduction or T-dualities along the half 
spacelike half timelike direction makes sense. Of course one might try to do the 
TbT operation along flat space-time directions but this will not generate any new 
solutions. In the literature TbT was generally done on asymptotically flat directions, 
see for example and references therein, however, when one tries to do the TbT 
operation along very curved directions as in our case, one will generally get into the 
same problem as we do. 

The above issues tell us that the duality chasing idea may be a futile endeavor 
to look for chemical potential from the gravity side. There exists a much better way 
to determine the source of chemical potential involving considerable less effort. We 
turn to this in the following section. 



3.3 Chemical potential from D7 brane world-volume gauge theory 

Another way to introduce chemical potential is to study the U{1) gauge field on the 
D7 brane as in P] (see also [^). To study this case thoroughly one has to take into 
account not only the backreaction of the D7 brane but also the backreaction of the 
gauge field on the brane. 

We first solve for the gauge field. To do this, we consider the solution from p 
where we already took into account the backreactions from the seven-brane as well 
as from the fluxes and black hole. This way the metric ansatze will be more generic 
than what we took here. Of course in 0, the warp factor and the hi coeffcients of 
the internal metric were left undetermined, so the chemical potential that we will get 
using this background will be left in terms of these coefficients. Our aim then would 
be to plug in the values of the coefficients that we determined here by considering the 
backreaction from the black hole and compute a representative chemical potential. 

Our starting point then is the metric^^: 



gttdt^ + dx\2z) + y/h{grrdr'^ 



r^dMl) 



(3.10) 



where dM^ is the same as ds\ in ( |A.1D . We take the simplest D7 brane embedding 
profile which is 02 = 0, ^2 = 6'2(r), so the induced metric for D7 brane is 

dsl = —={—gttdt^ + dxl2z) + ^ i.9rr + r'^hsh4{dr92)'^)dr'^ + r^hi^dip + cos^i(i0i)^ 
yh L 
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Note that in this section we are taking a metric ansatze more general than the one considered 



in (2^). Therefore we will assume that the radial range is r > r/j and not ( 2.31 ). Of course the 
metric of the internal space could still change to 0(r^/r"') but we will not address the issue here. 
A more elaborate analysis will be presented elsewhere. 
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+r'^h2{d6\ + sin^ Oid(j)l) + cosip{dr62)drd6i + smip sm6i[dr62)drc 



(3.11) 



In the UV region the NS and RR two-forms, B2 and C2 respectively, can be neglected. 
If we allow the seven- brane embedding to be wholly in Region 3 of [0, then the 
scenario will be no different from the ones in P|. Therefore by the same argument 
in we can only turn on the At component on the brane which will be a function 
of the radial coordinate. The DBI action for the D7 brane is, 



I = Td7 d®xe~Vv/M^sin^i 



gtt9rrh2 + gtth2dr(l)'dr(l)'Vh + /is^^^^ 



+gur'ih2hh^- hi/ A)idre2] 



1/2 



Neglecting the fluctuations in 0* i.e keeping dr4>^dr 
At is 



a 



y/hjl^ h2 Or At 



V9ugrrh2 + h2Fl + gttr^{h2hh^ - hl/A){drd2Y 



(3.12) 

0, the equation of motion for 
-- drD = (3.13) 



At r — )■ 00, one can find At = fi + ^ where /i is the chemical potential on the field 
theory side. And /i can be described as 



^•0 



°° , I gugrrh2 + gttr\h2h;h^ - /i^/4)(ag2)^ 
e-^^r^hihl-h2D^ 



(3.14) 



where tq » with tq to be a point in Region 3 determined by the minimum dip 
of the seven-brane in Region 3. For simplicity, we have also assumed the chemical 
potential for ^('"o) = 0. 

Note that the chemical potential in ( |3.14j ) depends on the precise embedding 
62{r). To eliminate 62{r) in ( ^.141 ) we need to determine the equation of motion for 



this. For this we first construct the Legendre transform of eq. ( 3.12 ) with respect to 
D to eliminate At in the action. The new action is, 



X = Td7 j d^xe-'t'r^ sin 9i 



hih2 {gugrrh2 + gttr'^{h2h^hi - hl/A){dr92Y) 

1 



(3.15) 



With this the equation of motion for 92{r) is, 

gtt{h2h^h^ - hl/A)dr92 



hih2 



D2 



\/gttgrrh2 + gttr'^{h2h^hA - hl/A){dr92) 



drK = 
(3.16) 
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After simplification we find that the chemical potential can be expressed in terms of 
two constants D and K in the following way: 



^i= I dr F^,{r-{hi}) (3.17) 

where we have defined 



ro 



(r; {h,}) 



(3.1^ 



The above formula for the chemical potential although elegant, is unfortunately not 
very illuminating for us because the warp factors hi are now all in Region 3 whereas 
in this paper we have concentrated ourselves only in Region 1. It could be that 
the hi continue to remain the same from IR to UV, but we don't have any a-priori 
information on this. In addition to that we expect the seven-brane to go all the 
way to r = r/j in Region 1. Therefore a more generic formula for the chemical 
potential, that takes into account various regimes of validity of our construction, can 
be expressed succinctly as^^: 

POO 

/i= / dr T^{r;{h,}) (3.19) 

= dr (r; {h^^}) + H dr 7,^^ (r; {hf}) + f dr 7,^^ (r; {hf}) 

where h^^^ are the warp factors of the internal space in Region n with the corre- 
sponding function [r] {/ij^^^jj. One interesting thing about (|3.19|) , and also of 
( ^.17| ), is that the warp factor h doesn't appear in the final formula (at least for the 
embedding that we have chosen). We also note that, in the limit: 



= 0, « log I ^ ^ ' I (3.20) 

which means that the value of the dilaton at any point in the radial direction is 
bounded above by the log function, the formula for the chemical potential drastically 
simplifies to: 

lj = D \ I dr —= + I dr —= + I dr —= (^.21' 



r^h'^^'JhY' •''•-i" r^h'^^^'JhY' "^'■0 r^h'i'^hX 



•^^In deriving this formula we have ignored the puU-backs from the NS and RR three-form fluxes 
from Regions 1 and 2. It will be interesting to see how the chemical potential depends on these 
informations. Here, however, we will suffice ourselves with approximate estimates, and a more 
detailed exposition will be addressed elsewhere. 
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The actual value of the chemical potential in this limit now depends on the value 
of the dilaton in the three regions as well on the internal warp factors To 
estimate the value of the chemical potential let us now consider a toy example where 
the dilaton is approximately constant over the three regions and the internal warp 
factors remain same in Regions 1, 2 and 3, i.e: 

^ gs, D = corl (3.22) 
1 9' 2 6 ^ 1 + ^ 6 ^ 



where cq is a constant and the scaling of D is motivated from f^. Plugging (|3.22|) in 
( p.21|) we get our final estimate for the chemical potential as: 



/i = - dcoThOs (3.23) 

where the sign of n will be determined by the sign of Cq. Note that our estimate says 
that the chemical potential is very small (provided \co\ ^ 1) and is proportional to 
Th, the horizon radius. In the dual field theory, the temperature is related to the 
horizon in the following way 

T = ^'^J^ (3.24) 
47r^h{rh) 

This means that in terms of temperature, the chemical potential is given by 

/i = coTrl^/h^ (3.25) 

where we have taken the appoximation that g{r) = 1 — r^/r'^- As seen from our 
numerical analysis, (^o ^ 1 and we can effectively ignore the modification of g{r) due 
to the presence of fluxes, while Aq near the horizon is not negligible. Thus we have 
to consider h{rh) = h^{rh) + h^{rh) in the above expression. From (|3.25| ), we see that 



chemical potential is approximately linear in terms of temperature, however as 




iVW^) - L\ll + %|^log(^) + Ao(r,) (3.26) 



with r^, being another relevant scale, ( p.25|) also involes terms ~ log(T). But this 



is precisely consistent with the logarithmic running of the gauge coupling with tem- 
perature and hence our analysis reveals that chemical potential is sensitive to the 
underlying structure of the gauge theory, which by construction is non-conformal. 

4. Discussions and Conclusions 

In this paper we have addressed one of the crucial issue left undetermined in our 
earlier papers T^, rH|, namely the backreactions from black hole, branes and 
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fluxes on the background geometry and on the various UV completions. We found 
that in certain hmit (|2.28|) the background EOMs aflow for a perturbative suppres- 
sions ( p.29p under which the backreactions from branes and fluxes can be consistently 
ignored. This means we only had to study the backreaction from black hole on the 
geometry, a task, which in the same limit, simplifles immensely to a set of three 
second-order partial differential equations ( |2.15D , ( [^.33| ) and ( [^.45D provided we are 
in the radial range ( |2.31| ). Solutions to these equations are provided in Appendix 
ly. Of course the analysis that we present here is only meant for the IR regime of the 
theory, but it would probably be possible to extend this to Region 3 where UV caps 
were introduced in P, [l^. Two challenges stiU remain: one, to study the equations 
at r = r/i in Region 1 and two, to study them in the intermediate buffering region 
i.e Region 2. In both cases the analysis may get very involved because for the flrst 
case one would now have to solve all the twenty internal Einstein's equations; and 
for the second case the (p, q) flve-brane sources and fluxes will further complicate the 
scenario. These details are left for future works. 

We took few other directions in this paper too, many of which are not restricted 



by the constraint (|2.31| ). We gave a detailed analysis of the backreactions on the 
fluxes from the black hole and the flavor seven-branes. These backreactions tend to 
make the three-form fluxes non-ISD. Interesingly the effect of the horizon radius 
on the fluxes is, to the order that we present here, implicit: the appearance is via the 
resolution parameter F (or a^). The fact that the resolution parameter would have 
dependence on the black hole radius was anticipated in |^ |T^. Here we conflrm that 
prediction. 

Another direction that we took here is related to the study of the running cou- 
plings. Our analysis predicts that the two couplings change under various backreac- 
tions from the dual gravity side, but their RG runnings remain the same as for the 
non-thermal case. We believe this should be the correct way to interpret thermal 
beta functions. An alternative interpretation where the thermal beta function is as- 
sumed to be different from the non-thermal case is probably not a very meaningful 
conclusion. 

Our flnal analysis, ignoring the detour that we took to study dipole deformed 
quarks, is the study of chemical potential. We showed how the chemical potential 
should be studied in a theory with an underlying RG flow. Our analysis yields, by 
ignoring certain subtleties, a result of the form ( |3.19| ). Under simplifying assum- 
tions this gives us ( [3.231) which turns out to be a reasonably good estimate. A full 
analysis however will require us to compute, among other things, the internal warp 
factors h^[^\ Furthermore, going beyond QsNj and ggM'^/N orders, and estimating 
the background would yield interesting results for the IR regime of large N thermal 
QCD. All of these are of course very challenging questions but, as we saw in this 
paper, there may exist corners of the solution spaces where one may flnd unexpected 
simpliflcations. It is worth at least, if not for anything else, to look for these regimes 
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of simplicity in a subject that is notorious for its unyielding complexity. 

Acknowledgement 

We would like to thank Miklos Gyulassy and Guy David Moore for helpful discus- 
sions. The work of M. M is supported in part by the Office of Nuclear Science of 
the US Department of Energy under grant No. DE-FG02-93ER40764, the work of 
K. D is supported in part by NSERC grant, the work of F.C and P. F is supported 
in parts by Schulich scholarships and FQRNT grants. S. V would like to thank the 
High Energy Group at McGill University for financial support. 



A. An exercise in duality chasing 

In this section our aim would be to understand how to switch on a gauge field, Aq, 
along a five- dimensional space in type IIB string theory. The five-dimensional space 
will be the dual gravitational description of a gauge theory with running couplings. 
This would mean that all the type IIB fields, namely the three and the five-forms 
including the axion-dilaton are all switched on. We will generate this solution using 
two stages of duality transformations. The first stage of duality transformations 
would convert a simple torsional background to another background that has both 
the three-forms including five-forms and axion-dilaton switched on. The second stage 
of duality transformations would further convert this background to the one that we 
want. 

Our starting point would then be to take the following torsional background: 

dsjo = hh^hodsl + h-^e^dsl (A.l) 
where is the dilaton and we have defined ds^ and ds\ in the following way: 
ds ^ — — cjiidt -\- dx-^2^ 

dsl = hidr'^ + h2{d'ip + cos9id(f)i + cos62d(f)2f' + h^idOl + sin^ Oidcpl) + h^^iJiQ d6\ + sin^ d2d(j)l) 
+hr, cos ip {d9id92 — sin 9i sin ^2(^0i02) + ^5 sin ip (sin 9id92d(f)i — sin 92d9id(f)2) (A. 2) 

with (/i, hi) being the required warp factors, and the torsion H-^ is given by the usual 
relation: 

= * d{e-^*J) (A.3) 

J is the fundamental form that we will compute separately in Appendix Now 
using our first stage of duality transformations, as described in details in , we can 
easily get the following background from ([A.l|) , ( [A.3|) and the dilaton 0: 



F^ = hhl cosh Pe^'f' * d(e-2'^ J) , 
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= -hhlsmhf3e^'^d{e-^'^J), 
F5 = + i^)dAo A A dx'^ A dx'^ A dx^, 

ds^ = hodsl + hidr"^ + h2{dip + 008 6*1^01 + 0086*2^02)^ + h^^dOl + sin^ Oidcpl) + h^^iJiQ dO^ 

+ 8in^ 02d(pl) + /15 C08 -0 {d9id02 — 8in 6*1 8in 6*2(^0102) + /is 8in -0 (sin 9id92d(pi — sin 02dOid(j)2) 

(A.4) 

For this solution to be supersymmetric, the minimum requirement, in the absence 
of a black hole, is hh^e^'^ = 1. Using this we can compute the NS i?-field from the 
three-form = dB2 as: 

B2 = -e-^'^sinh/3J 

= -e"^"^ sinh l3{^/F\F^dr A {d^p + cos6'id0i + cos 612^02) 

sin -0 sin ^1 sin^2('U^ + za)d(j)i A d(j)2 ~ ABzsm.6id(j)i A dOi 
+AB cos -0 sin 62(10 + za)d9i A (i02 

-[.450(1(7 + za) - C^F^ - a^A^{z - wd)] sin02C?^2 A (i02 

-[ABaz - C\/ F^- a^A^w\ cos'0sin6'irf^2 A dc^x 

+ [ABaz - C^F^ - a^A^w\ s\ni\}de2 A dQ^) (A.5) 

where A, 5, C, a, w, J are worked out in Appendix p. 

In our second stage of duality on ( |A.4D , we will do a T-duality followed by a Boost 
and then another T-duality (TbT) along 02 direction. To make the TbT operation 
simpler, we can do a gauge transformation to eliminate the B^-^i components in the 
following way: 

B2 = B2 + dCu with Ci = j B^,id(t)2 (A.6) 

Now that the B2 field has no B^^i components, TbT is straightforward. The resulting 
metric is: 

ds^ = ( — — _ — _ — ^^2 _|_ dxf^^ + hidr'^ 



A cosh (3 — hoQu sinh (3 



Acosh^ (3 



+ \ h2-hl COS^ 92A{1 -y- ; -2—) (# + COS 9id(j)i + COS 6'2rf02)^ 

V A cosh p — h^gu smh p 

+h3{del + sin^ 6*1^01) + h^ih^ dOl + sin^ 6*2(^02) 

2 



A cosh^ /3 
^ cosh^ 13 — h^Qu sinh"' (3 



( 1 ^ ^-"2 ^ — 2~^ ) {h2h^sm.'ip cos 2 sm. 6 2d'iljd6 1 



hi siv? -0 sin^ 62 ,„2 ^5 "^o^^ ip siv? 61 siv? 62 — h2h^ cos ip sin 9i cos 61 sin 6*2 cos 62 , , 2 
dui — "01 



— {h\ sin'' 6*2 + 2/12/14 sin^ 6*2 cos^ 62)d(t)\ — (2/2,2 cos cos 62 — cos -0 sin 9i sin 6*2) 
X — — — -d(f)id6i + /;-2/i5 cos -0 sin 61 sin 6*2 cos 62d'ipd(f)i + /ij cos 6*1 cos^ ^2*^0 
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A cosh (3 

A cosh^ 13 — hogtt sinh"' (3 
A cosh (3 



— f 1 K K — I (h2 cos6'i cos 92d(bi(b2 + /i2^4 cos6'2 sin^ 92Adihd(h2) 



-5 2 — \ ^5 '^os ^1 92d(f)2d(j)i + /15 sin sin 92d(p2d9i 

A cosh /3 — /loS'tt sinh /3 

cosh. y5 

H 2 2 — ^2 "^os^ ^2(1 ~ cosh [3)dipd(l)2 + /i5(cos ipd9id92 + sin sin 9id92d(f)i) 

A cosh ^ — /io5'tt sinh /9 

(A.7) 

where (3 is the boost parameter, and we have defined A as: 

The metric (|A.7| ) is accompanied by all the type IIB RR and NS form fields. These 
can be easily worked out. Here we will just quote the B2 field and the dilaton (p: 

= 0- ^log(cosh^ (3 - hoQu sinh^ /3) + hog A 

{A - hogtt) sinh P cosh P n ■ u n 

Bq^^ = -y- ; — t:; , = Ah2 COS 92 smh (3 

A cosh (3 — hogtt smh (3 

A A 

Bo^i = — (2/i2 cos 9i cos 92 — h^ cos ip sin 9i sin 6^2) sinh /3, Boe^ = — h^ sin ip sin ^2 sinh [3 



(A.9) 



It is now easy to see what the five- dimensional gauge-field would look like. Its simply 
the following wedge product: 

Ao = Bo^^d(j)2 + Bo^dip + Bo4>^d(j)i + Boeid9i (A. 10) 

Notice that for the metric and dilaton to be regular we require the metric component 
gtt and the warp factors to satisfy: 

hogtt{h2 cos^ + /i4 sin^ 92) ^ coth^ /3 (A.ll) 

If we demand this requirement to hold for all the values of ^2 and (3, then we need 
the following two scenarios: 

If hi > h2, then hoh^gtt ^ 1 

If h2 > /i4, then /io/i2fi'tt ^ 1 (A. 12) 
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B. The fundamental form J 



We first define the vielbeins as follows, 

^ — s/h^idi^ + cos 9id(f)i + cos 92d(f)2) 

— A{sim/jsm9id(f)i + cos'^dOi — ad92) 

— b(^w sin ^2C^02 + -z(cos ip sin — sin ipdOi + a sin ^2O?02) 
= Clzsm92d(p2 — w{cos'ijj sm9id(f)i — simjjd9i + a sin 6*2(^02) ) (B-l) 



where A, B, C, a, w, z must satisfy, 

h 3 

B''z^ + CW = h^, B''z-C^w = hs, B^w^ + C^z^ = h^-^ (B.2) 

The fundamental form J is, 

J = A + A + A 

= -s/hih^dr A (d'^ + cos + cos ^2C^02) 

sin sin 9i sin ^2('Ji' + za)d(j)i A d02 — sin 9id(j)i A d^i 
cos t/j sin ^2('Ji' + za)d9i A d02 

— [74Sa(t(; + za) — C-s/hJiQ — a?A^{z — wa)\ sin ^2C^^2 A d02 

— [ABaz — C-\/hJiQ — aP'A'^w] cosipsm9id92 A d4>i 

+[ABaz - C^/hJiQ - a?A^w] sin ^'^^2 A d9i (B.3) 
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C. More details on the squashing and the warp factor com- 
putations 

In the absence of blackhole that is e^^ = 1, the warp factor a = e^^ = 1/h only 
depends on r, 6i and 62 even when D7 back reaction is taken into account by con- 
sidering the running axion-dilaton field In this extremal limit Vh = 0, we have 



ISD three-form fluxes G3 and the internal metric gmn describes a Ricci flat deformed 
cone. For the non-extremal case, we will demand similar behavior for the warp fac- 
tor h and will find that such solutions do exist. Using h = h{r,6i,62) only, in the 
non-extremal case we get, 

dh = d(^[drhC + e-'^ {f^'^de.h V, + f^'^de^h r/^)] ^^^^^^^^^sin^isin^^) 
= dV (C.l) 



where we have used our metric ansatz (|2 . 1| , |2.3| ) and definition of the five-form flux 
( p.43|) . We have only kept linear terms in F,Q and this is justified as we look for 



solutions F, ^ ^ 1 and ignore terms higher order term. In the above we have also 
defined 



9pq = e-'%q (C.2) 



where p, q run over the compact directions and thus ( p.2|) implies gpq is indepndent 
of B. Here C,''7i are five- forms given by 

C = dip A dcpi A (i02 A dOi A c/6'2 
rji = dip A dr A dcpi A d(p2 A ^6*2 

r]2 = dip Adr A dcpi A d(p2 A dOi (C.3) 

Just like the extremal case, we will assume that de.h ~ 0{glNfM'^/N) and we 
will find that this choice is consistent with all the Einstein equations and equations 
for the fluxes . With this assumption, we readily get upto 0{gsM'^ /N), and ignoring 
0{F,G,g)0{gsM^/N) (since F, G, < 1) 

f'''de,h = -g'!;'^de,h (C.4) 

where i = 1, 2 and g^'^ is zeroth order in M, Nf. But at zeroth order in M, Nj, the 



compact five dimensional internal space A^s is exactly the deformed cone and thus 
^q'' is precisely the metric of deformed T^'^. Our ansatz for the black hole factor e^^ 
is given in (|2]4|) where G is at least 0{M/N, g.M"^ /N, g^Nf). This is a sufficient 
condition as in the absence of five-branes and seven branes, we have AdS x T^'^ 
with black hole where e"^^ = 1 1 and fh ^ b. This is because for large r, the 
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deformed cone becomes the regular cone and considering fh 3> a, we are effectively 
putting a black hole in a regular cone. In other words, the non-extremal limit of the 
geometry only 'sees' the regular cone and deformation of the cone is hidden behind 
the black hole horizon. This also means our non-extremal solution is valid only for 
large horizon, that is the non-extremal solution only captures the large temperature 
deconfined chirally symmetric phase of the gauge theory. The extremal solution 
without any black hole is dual to the confined phase. 
Now using (|]4|) in V reads 



V 



1 



drhC + -rigrde^hr],) 



' h 

7^ 



+ 



1 - 

[drh\ + g'/^deM' Vi\ 



^. sinMng, (l + F + G/2) 



108 



• sin6'isin^2 



1 - ^ 



108 

sin6'isin6'2 
108 



drh^ C + 



1 — 



(F + C;/2)smMn^2 

108 ^ 



(C.5) 



where we have only considered upto 0{gsM'^ /N) terms. Here is the Ouyang 
solution and is the correction due to the black hole which alters the internal 
compact space M.^. But the Ouyang solution satisfies Bianchi identity exactly as: 



d 



{drh\ + g'/^de,h\) T 



jSin^^isin6'2 
108 



(C.6) 



Using this in (|2.37|) we get 



drh\+--g'/^de^h\+'^^ 
9 9 



sin^isin^a 5 (F + G'/2) sin^isin^a 



108 



108 



drh\ 



which gives us ( |2.45| ). The derivations of ( ^.151 ), ( [^.33| ) and ( [^.35| ) have already been 
discussed in section 2.2. 

We will now solve the four equations (p.45|) , ( p. 151) , ( p.33| ) and ( p. 351) by ignoring 
all terms of 0{gsNf). In this limit, all angular dependences vanish and all the func- 
tions A, F and G are only functions of the radial coordinate r. This also means we 
are ignoring the back reaction of the seven branes and our solution should be consid- 
ered as the non-extremal generalization of Klebanov-Strassler theory with modified 
UV behavior. For Nf = 0, with e~'^^ = h = + h^, v/e take the following ansatz 

= — {Mr) + ^i(r)log r + A2(r)logV) 



.25 



4 + 9o{r) + fi'i(r)log r + 5(2(r)logV 



F = Fo(r) + Fi(r)log r + F2(r)log^r 
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(C.7) 

With our ansatz, only taking upto linear order terms in Ai, Fi and gi one obtains that 
the equation derived from (|2.35|) is trivial. Also upto linear order, Ai = A2 = Fi = 
-^2 = S'l = ^72 = is a solution with Aq.Fq.Qq being the only non-trivial functions. 
The equations resulting from ( |2.45| ),( p.l5| ) and ( p.33| ) are as follows 



ii) rA'o' - 3A'o - = 

{lii) — fl + 56r^go + 16r^Fo + 4rf^A'o + idr^g'^ + 24r^F^ + 12rr^F^ 
Ntc 

+ 7r'g'^ + Ar'F;; - ArhlF^; = ( 



To solve these second order differential equations, all we need to do now is 
specify the boundary conditions. As we have second order differential equations, we 
can choose two boundary conditions. A priori we do not know where the horizon is, 
that is we do not such that e'^^'^'^h)=o ^ cannot specify the boundary condition 

at the horizon. Additionally we cannot take r to be smaller than the range ( |2.31| ). 



However, since we are looking for solution such that asymptotically we recover the 
extremal geometry, we can impose the following boundary conditions: 

lim^^oo AQ{r) = 
lim^^^oo goir) = 

lim^^oo Foir) = (C.9) 



From the form of the equations in (|C.8|) , we see that inverse power series in r is a 



possible candidate for the solutions that obey the boundary conditions ( |C.9| ). On 
the other hand, as already discussed in section 2.1, we expect Ai,gi and F^ to be 
proportional to the horizon Vh- Thus our anstaz is 

Cl(-f)\ (CIO) 



where a^, and Cfc are atleast 0{M/N,gsM'^/N), and the radial coordinate r is 
assumed in the range ( p.31| ). The boundary condition (|C.9|) implies 

4 = fo= Co = (C.ll) 

We can further choose three other boundary conditions. Again since (a) we do not 
know where the horizon is and (b) the radial coordinate is constrained by (|2.31 ), we 
will choose the following boundary conditions: at r = 00 and choose 

lim,.^oo A^'oir) = 



- 51 - 



limr_>oo g'oir) = 

lim^^oo i^oW = (C.12) 

which is automatically solved by our ansatz ( |C.10|) . With the set of boundary con- 
ditions ( |C.9| ) and ( |C12| ), we solve ( |C.8| ) numerically. The exact solution (whose 
validity should be considered for r > (N/MY^^rh) is plotted in Figures 1, 2, 3. 
Observe that the numerical solutions are consistent with the analytic behavior in 
( P-lOp . When M = 0, equations ( |C.8|) imply that we have the trivial solution, i.e. 
^0 = go = Fq = 0. But since M 7^ 0, we must have non-trivial solutions to satisfy 
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